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SUMMARY 


The  initial  phade  of  this  work  was  a feasibility  study  of  an 
- HCl  mixing  gas  dynamic  laser.  The  results  of  that  study  indicated  that 
a high  energy  laser  system  based  on  the  - HCl  system  was  Indeed  feasible. 


This  laser  device  concept  envisions  the  use  of  a aolid  propellant  to  generate 


high  temperature  H'J  - HCl  mixtures  which  are  expanded  supersonically  to 


create  the  proper  conditions  for  efficient  H2  ^HCl  VQ»-  V transfer  and  lasing 
on  HCl.  This  type  of  laser  device  has  several  significant  potential  advan- 
tages as  a high  energy  device  for  military  application,  l.e.  , the  use  of 
compact,  field  storable  solid  fuel  technology,  high  specific  power  (80  - 100 
Kj/lb),  and  lasing  of  HCl  In  the  A.CJ^^a^ospherlc  window.  >'  ' 

The  objective  of  the  current  phase  was  to  experimentally  verify 
the  V - V and  V - T kinetics  for  the  important  H2  - HCl  processes  and  to 
develop  solid  propellant  gas  generators  that  could  be  used  to  generate  the 
required  gases  at  attractive  weight  yields.  The  results  of  the  current  phase 
have  verified  that  the  kinetics  are  capable  with  operation  of  a H2  - HCl 
laser.  In  addition,  the  gas  generator  development  program  has  succeeded  In 
producing  the  required  gases  at  weight  yields  that  make  the  system  attractive 
for  Array  field  application. potential  of  even  higher  gas  yields  that 
would  make  the  system  even  more  attractive  Is  also  possible  but  requires 
additional  propellant  development.  The  fundamental  remaining  problem  area 
Is  the  effect  of  H- atoms  on  performance. 
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1.0 


INTRODUCTION 


The  initial  phase  of  this  work  was  a feasibility  study  of  an  • 
H2-HCI  mixing  gas  dynamic  laser.  The  results  of  that  study  indicated  that 
a high  energy  laser  based  on  the  H2-HCI  molecular  system  was  indeed  feasible. 
This  laser  device  concept  envisions  the  use  of  a solid  propellant  to  generate 
high  temperature  H2-HCI  mixtures  which  are  expanded  supersonically  to  create 
the  proper  conditions  for  efficient  H^-^HCl  V-+-V  transfer  and  lasing  on  HCl. 
This  type  of  laser  device  has  several  significant  potential  advantages  as  a 
high  energy  device  for  military  application,  i.e.,  the  use  of  compact,  field 
storable  solid  fuel  technology,  high  specific  power  (80-100  Kj/lb),  and 
lasing  of  HCl  in  the  4. Op  atmospheric  window. 

In  the  initial  study  phase  the  critical  kinetic  and  propellant 
issues  that  required  further  research  were  identified.  In  the  kinetics  area, 
it  was  determined  that  the  three  most  critical  kinetic  rates  that  control  the 
laser  performance  were  (i)  the  H2/HCI  V->V  process  for  which  no  data  were  then 
available,  (ii)  the  HCI/H2  V->T  process  for  which  only  a room  temperature  data 
point  had  been  determined,  and  (iii)  the  HCl/H  V-+-T  process  for  which  only 
room  temperature  measurements  existed.  In  the  propellant  area,  it  was  deter- 
mined that  gas  generators  with  the  desired  gaseous  weight  yields  and  gas 
composition  must  be  developed. 

The  purpose  of  the  present  investigation  (Phase  II)  was  to  provide 
experimental  data  on  the  rates  of  the  H2/HCI  V->-V  and  HCI/H2  V->^T  processes  at 
temperatures  in  the  range  of  1000-2500°K,  to  develop  the  gas  generators  with 
the  required  weight  yield  and  gas  composition  characteristics  as  well  as  to 
perform  additional  analytical  studies  for  the  design  of  an  appropriate  proof 
experiment.  The  results  of  this  work  are  reported  in  this  docviment. 
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2.0 


KINETICS 


2 . 1 Background 

As  was  emphasized  in  the  early  stages  of  the  H„-HC1  GDL  feasibility 

1 ^ 
study  , the  vibrational  energy  transfer  kinetics  play  a critical  role  -in 

determining  potential  laser  performance.  The  vibrational  temperature,  gain, 

and  extractable  power  are  all  highly  dependent  on  the  rates  for  V-T  and  V-V 

exchange.  The  principal  vibrational  energy  transfer  processes  used  in  the 

harmonic  oscillator  model  are  listed  in  Table  2.1. 

There  is  presently  a large  body  of  experimental  data  on  process 

2-5 

(1)  obtained  from  several  different  experiments  and  covering  a large  tem- 
perature range  from  150  ^ T _<  2100.  The  rate  constant  for  process  (2)  has 
been  measured  at  room  temperature^ , and,  quite  recently,  some  shock  tube- 
laser  fluorescence  data  by  Bott  have  become  available  which  provide  some 
additional  information  about  process  (2)  up  to  T = 770“K.  These  most  recent 
data  by  Bott  will  be  discussed  later  in  this  section. 

Experimental  data  for  V-T  process  (3)  exist  from  three  separate  ex- 

g 

periments  - a low  temperature  experiment  covering  the  temperature  range  ftom 

9 10 

approximately  50 ®K  to  450°K  and  two  high  temperature  shock  tube  experiments  * 
covering  the  temperature  range  from  1100°K  to  2700®K.  These  data 
allow  us  to  make  a reasonably  good  curve  fit  to  the  rate  constant  for  process 
(3)  over  the  temperature  range  of  interest  for  the  proposed  H2-HCI  GDL. 

While  there  is  at  present  no  experimental  data  in  the  literature  for 
process  (4)  , it  is  not  expected  that  this  reaction  will  make  a significant 
contribution  to  the  vibrational  energy  relaxation.  According  to  theories  of 
the  SSH  type,  HCl  is  expected  to  be  less  efficient  at  deactivating  vibrationally 
excited  H2  than  H2  itself.  This  fact  coupled  with  the  small  concentration  of 
HCl  relative  to  H2  in  the  laser  medium  make  the  contribution  of  process  (4) 
negligible. 

Perhaps  the  most  critical,  yet  most  poorly  determined  rate  is  that 
for  the  V-V  exchange-process  (5).  Previous  to  the  present  experiment,  no  exper- 
imental data  existed  for  the  rate  of  process  (5)  and,  therefore,  the  laser 
modeling  calculations  have  had  to  be  done  using  an  estimated  rate  expression 
for  this  process  derived  from  a combination  of  empirical  correlation  and  theory^ 
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TABLE  2.1 


H^/HCl  VIBRATIONAL  ENERGY  TRANSFER  PROCESSES 
HARMONIC  OSCILLATOR  MODEL 


V T Processes 

HCl  (v  = 1)  + HCl  HCl  (v  = 0)  + HCl  + 2886  cm^  (1) 

HCl  (v  = 1)  4-  HCl  (v  = 0)  + H2  + 2886  cm"^  (2) 

H^  (v  = 1)  + H^  H2  (v  = 0)  + H2  + 4159  cm~^  (3) 

H2  (v  = 1)  + HCl  H2  (v  = 0)  + HCl  + 4159  cm~^  (4) 

V ->  V Process 

H2  (v  = 1)  + HCl  (v  = 0)  ^ H2  (v  = 0)  + HCl  (v  = 1)  + 1273  cm“^  (5) 
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On  the  basis  of  the  above  analysis,  it  is  clear  that  the  Important 
molecular  vibrational  rate  constants  that  required  experimental  investigation 
are  those  for  the  V-V  process  (5)  and  the  V-T  process  (2) . The  present  shock 
tube  experiment  was  designed  to  make  measurements  of  these  rate  processes 
over  the  temperature  range  of  roughly  1000*K  to  2500“K. 

2.2  Experimental  Technique 

The  experiments  were  carried  out  on  a conventional  1.5-lnch  diameter 
shock  tube.  A photograph  of  the  general  shock  tube  facility  is  shown  in 
Figure  2.0. 

The  basic  technique  employed  monitoring  the  infrared  fluorescence 
emitted  by  the  HCl  behind  Incident  shock  waves.  It  has  been  shown  that  the 
Infrared  radiation  intensity  from  the  fundamental  vibrational  rotation  band  is 
proportional  to  the  vibrational  energy  in  that  mode.  Therefore,  following 
the  time  history  of  emission  from  the  fundamental  band  is  equivalent  to 
measuring  the  rate  of  change  of  the  vibrational  energy. 

The  shock  tube  was  constructed  of  1.5-inch  i.d.  stainless  steel  with 
a 4-foot  long  driver  section  and  a 15-foot  driven  section.  Shock  waves  were 
produced  by  bursting  scored  aluminum  and  cold  rolled  steel  diaphragms  with 
helium.  The  shock  speeds  were  measured  with  a series  of  six  thin  film, 
platinum  heat  transfer  gauges  mounted  at  various  positions  down  the  tube. 

The  output  from  the  heat  transfer  gauges  were  displayed  on  a speed  raster 
sweep  of  a suitably  modified  Tektronix  oscilloscope.  A typical  photograph 
of  such  a speed  raster  sweep  with  the  differentiated  outputs  of  the  various 
heat  transfer  gauges  is  shown  in  Figure  2.1e, 

The  infrared  emission  intensity  was  detected  behind  the  incident 
shock  by  means  of  the  optical  system  shown  in  Figure  2.2.  The  radiation  first 
passed  through  a 0.50-lnch  diameter  CaF2  window  contoured  to  match  the  inner 
radius  of  the  shock  tube.  The  light  was  then  collected  by  a 10-cm  diameter 
aluminum  mirror  with  a 30-cm  radius  of  curvature,  directed  to  a flat  aluminum 
mirror,  and  finally  Imaged  with  1:1  magnification  on  a 1-mm  diameter  detector 
element.  The  detector  was  a liquid  nitrogen  cooled  InSb  detector  operated  in 
the  photovoltaic  mode.  An  Infrared  filter  with  a bandpass  of  3.2  to  5.8v  was 
placed  before  the  detector  to  Isolate  the  fundamental  HCl  emission. 
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IR  INTENSITY 
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Figure  2.1  Sample  Fluores 


cence  and  Shock  Speed  Oscillograms 
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The  test  section  of  the  shock  tube  was  constructed  of  Lucite  to 
minimize  the  stray  light  that  might  reach  the  detector  by  internal  reflections 
off  the  metal  shock  tube  walls.  Before  this  Lucite  section  was  utilized, 
light  emission  was  detected  before  passage  of  the  shock  front  by  the  optical 
station,  presumably  from  radiation  emitted  down  the  tube  ahead  of  the  shock 
and  reflected  into  the  optical  system  by  the  back  wall  of  the  original  metal 
test  section.  Changing  the  test  section  to  Lucite,  which  is  a poor  reflector 
in  the  infrared,  eliminated  this  extraneous  light. 

The  output  of  the  IR  detector  was  amplified  by  a wide  band  low  noise 
preamplifier  of  20  Db  gain  before  recording  by  a 5A5B  Tektronix  oscilloscope. 

A 3 MHz  high  frequency  cut-off  filter  was  employed  Immediately  before  the  in- 
put to  the  final  amplifying  stage  of  the  oscilloscope.  The  time  response  of 
the  optical  detection  system  and  associated  electronics  was  approximately  Ip 
sec.  A sample  of  typical  IR  fluorescence  traces  is  shown  in  Figures  2.1a  and 
2.1b. 

Initial  fill  pressures,  p^^,  in  the  driven  section  were  measured  with 
a Wallace  & Tiernan  absolute  pressure  gauge  (0  - 200  torr)  and  an  oil  mano- 
meter. Typical  fill  pressures  ranged  from  about  3 torr  to  20  torr.  The  gases 
used  in  this  tudy  were  Argon:  Matheson  Purity,  minimum  purity  99.995%;  HCl: 
Linde  Specialty  Gases  L.E.D.  grade,  minimum  purity  99.995%;  and  H2*.  Linde 
Specialty  Gases,  ultra  high  purity,  minimum  purity  of  99.999%  and  a maximxim 
moisture  content  of  3 ppm.  The  argon  was  further  purified  to  remove  any 
possible  water  content  by  passage  through  a copper  coil  immersed  in  a dry  ice- 
acetone  bath.  The  HCl  was  subjected  to  a freeze-pump- thaw  cycle  between  77°K 
and  300®K  to  remove  any  residual  H2  that  might  have  formed  due  to  decomposition. 
This  HCl  was  then  passed  through  a copper  coll  at  dry  ice  temperature  to  re- 
move residual  water  vapor.  Finally,  the  H2  was  passed  through  a copper  coll 
immersed  in  a liquid  nitrogen  bath  to  trap  out  any  water  vapor.  The  gas  mix- 
tures were  prepared  in  a 4 liter  stainless  steel  mixing  tank  at  total  pressures 
up  to  100  psig  and  allowed  to  mix  by  diffusion  for  at  least  twenty-four  hours 
before  utilizing  them  in  the  experiments.  The  shock  tube  was  capable  of  being 
pumped  to  2 X 10  ^ torr  before  each  run  with  a combined  leak/outgas  rate  of 
less  than  7 x 10  torr /minute. 
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2.3 


Experimental  Results 


Since  the  rates  for  vibrational  relaxation  of  HCl  in  argon  are  well 
12  13 

established  " over  a substantial  temperature  range,  the  first  data  were  taken 

on  HCl-Ar  mixtures  for  the  purpose  of  calibrating  our  system.  As  expected, 

the  time  resolved  HCl  fluorescence  behind  the  incident  shock  was  found  to 

follow  a single  exponential  behavior.  Figure  2.3  is  a semi -logarithmic  plot 

of  the  normalized  HCl  emission  intensity  versus  lab  time  for  a 10%  HCl-90% 

Ar  mixture,  an  Initial  pressure  7.5  torr,  and  a shock  speed  of  1.17  mm/ysec. 

Fig.  2.4  summarizes  the  results  obtained  for  a 5%  HCl-95%  Ar  mixture.  The 

data  have  been  reduced  to  particle  relaxation  times  referenced  to  a post 

shock  pressure  of  1 atmosphere  for  various  post  shock  temperatures.  The 

large  Ar  dilution  assures  an  essentially  constant  post  shock  temperature 

during  the  relaxation.  The  solid  line  represents  the  results  obtained  by 
13 

Seery  for  a similar  mixture,  and,  as  can  be  seen  there  is  good  agreement 
between  our  data  and  that  of  Seery.  Data  obtained  in  other  HCl-Ar gon  mixtures 
demonstrated  similar  good  agreement. 

With  the  HCl-Ar  runs  completed  and  confidence  in  our  diagnostic 
system  and  gas  handling  procedures  established,  we  proceeded  to  make  shock 
tube  runs  for  various  HCl-H2-Ar  mixtures . 

The  vibrational  energy  relaxation  equations  for  a binary  gas  mix- 
ture of  harmonic  oscillators  has  been  derived  previously^'^  and  will  be  utilized 
here  without  further  demonstration.  The  specific  equations  used  are  as 
follows : 


dt 


I 'P/t 


(6) 


- V 

1 - exp(-e^) 


1 - exp(-6g) 


[1  - exp(-e^)] 


and 
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Figure  2.3  HCl-Ar  Shock  Tube  Run.  Semilogarlthmic  Plot  of  HCl 
Vibrational  EnerEV  Relaxation  Behind  Shock. 
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where  X = (E  - E)/E,  E Is  the  equilibrium  vibrational  energy  at  translational 
temperature  T,  and  E is  the  vibrational  energy  at  the  local  vibrational  temper- 
ature T^;  Ip  and  1 - i|)  are  the  mole  fractions  of  species  A and  B respectively; 
and  0 = hcm/kT.  Denoting  H2  as  species  A and  HCl  as  species  B,  x^,  x^g, 

Xba>  ‘'^ab  the  characteristic  time  constants  for  process  (3),  (4),  (1), 

(2),  and  (5),  respectively  of  Table  2.1. 

While  the  above  equations  may  be  linearized  to  provide  approximate 
solutions^^’  , for  the  present  analysis  solutions  were  generated 
numerically  by  computer.  With  Input  of  Initial  fill  pressure,  temperature 
(295®K),  gas  composition,  and  incident  shock  velocity,  the  computer  program 
calculates  the  post-shock  conditions  and,  depending  on  the  values  chosen  for 
the  unknown  rate  constants,  predicts  normalized  temporal  profiles  of  both  the 
HCl  and  H2  vibrational  energies  following  passage  of  the  incident  shock. 

From  a number  of  initial  computer  simulation  runs  of  various  exper- 
imental conditions,  it  became  apparent  that  determination  of  the  V-V  rate  for 
process  (5)  from  the  experimentally  observed  relaxation  profiles  would  require 
a reasonably  accurate  prior  knowledge  of  the  rate  for  V-T  process  (2).  Further, 
it  was  determined  that  there  was  a regime  of  experimental  conditions  where  the 
rate  for  process  (2)  could  be  extracted  essentially  Independent  of  an  evalua- 
tion of  the  V-V  rate  constant,  namely,  at  relatively  low  H2/HCI  concentration  ratios, 
i.e*.  , 1:1,  2:1  and  3:1,  and  at  the  lower  end  of  the  experimental  temperature 
range.  Under  such  conditions,  energy  sharing  between  HCl  and  H2  is  minimal  and  the 
relaxation  profile  is  only  very  weakly  dependent  on  the  V-V  rate.  This  fact  is 
illustrated  by  Figure  2.5  which  is  a calculation  of  the  temporal  profile  of  the 
HCl  vibrational  energy  behind  the  incident  shock  for  a 10%  HCl-20%  H2-70%  Ar 
mixture  at  an  initial  pressure  of  10  torr,  an  incident  shock  speed  of  u^  = 

1.0  mm/ysec,  and  a post  shock  temperature  of  915”K.  Comparison  of  the  two 
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Figure  2.5.  Shock  Tube  Calculation  of  HCI  Vibrational  Energy  Relaxa- 
tion. Shows  Regime  of ■ Temperature  and  H /HCI  Concentra- 
tion Ratio  Where  Relaxation  is  not  Sensitive  to  V-V  Rate 


calculated  curves  shows  only  a small  difference  in  the  predicted  relaxation 
behavior  even  though  the  assumed  V-V  rate  has  been  varied  a factor  of  100 
around  the  baseline  theoretical  rate. 

Such  computer  stimulations  made  clear  the  logical  sequence  of  ex- 
perimehts  to  be  performed  to  determine  both  the  HCI-H2  V-T  rate  and  the  HCI-H2 
V-V  rate.  First,  runs  would  be  made  for  mixtures  with  relatively  low  H2/HCI 
concentration  ratios  in  order  to  ascertain  the  rate  for  V-T  process  (2) . 
Following  an  accurate  determination  of  this  rate,  runs  would  be  made  for 
mixtures  of  high  H2/HCI  concentration,  where  the  relaxation  is  expected  to 
be  strongly  dependent  on  the  rate  for  the  V-V  coupling  process. 

A series  of  data  was  taken  on  a mixture  of  8.3%  HCl-16.7%  H2^-75%  Ar 
over  a temperature  range  from  1000“K  - 2000“K.  Figure  2.1a  is  a sample 
oscillogram  of  one  such  run.  A few  of  these  runs  were  analyzed  for  the  rate 
of  V-T  process  (2) , and  the  results  are  shown  as  the  square  data  points  in 
Figure  2.6.  A curve  fit  to  these  results  combined  with  the  established  room 
temperature  measurement  is  best  described  by  the  expression: 

k2  = 2.5  X 10  exp (-55. 9 /T^^^)  cm^ /particle”^  -sec  ^ (8) 

The  circles  shown  on  Figure  2.6  represent  the  data  of  Bott  and  Cohen^  between 
300“K  and  800°K.  On  comparing  these  two  independent  results,  there  appears 
to  be  a discrepancy.  Better  agreement  between  the  two  sets  of  data  is  possible 
if  one  takes  into  account  the  possible  contribution  of  the  H2-HCI  V-V  rate 
process  (5)  to  Bott  and  Cohen's  observed  HCl  relaxation.  Such  an  anlysls  of 
Bott  and  Cohen's  high  temperature  data  is  shown  in  Figure  2.7  and  described 
later.  As  can  be  seen,  assuming  a combination  of  a V-V  rate  based  on  our  theo- 
retical estimate  and  an  HCI-H2  V-T  rate  based  on  the  present  results  extrapo- 
lated to  770 “K,  one  can  predict  a relaxation  profile  that  follows  Bott's  ob- 
served decay  quite  closely  over  nearly  two  e- folding  times.  While  such  an 
interpretation  of  Bott's  data  is  premature  pending  the  experimental  determination 
of  the  V-V  rate  constant, it  does  provide  a preliminary  explanation  for  the 
apparent  discrepancy  between  the  two  experiments. 
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Figure  2.7.  Interpretation  of  Bott's  Observed  HCI/H2  Data  In 
Terms  of  V-V  and  V-T  Processes 
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Using  the  value  for  the  HCI-H2  V-T  rate  constant  established  from 
the  present  experiment,  runs  made  in  mixtures  of  high  H2/HCI  concentration 
were  analyzed  to  obtain  Information  on  the  value  of  the  V-V  rate  constant. 
Plotted  in  Figure  2.8  is  the  intensity  data  from  the  oscillogram  shown  in 
Figure  2.1b  taken  in  a 3%  HCl-57%  H2~40%  Ar  mixture  at  a post  shock  tempera- 
ture of  1170°K.  The  solid  lines  are  computed  curves  for  various  values  of 
the  assumed  V-V  rate,  the  baseline  value  being  the  theoretically  estimated 
rate  expression.  The  dashed  line  indicates  the  relaxation  predicted  if  the 
contribution  from  the  V-V  process  were  negligible.  It  is  clear  that  the 
data  plotted  in  Figure  2.8  indicate  a V-V  rate  constant  between  1/5  to  5 
times  the  theoretical  baseline  value  at  1170°K,  with  the  best  fit  shown  being 
that  utilizing  the  baseline  value.  Figure  2.9  is  another  similar  plot  for  a 
2.6%  HCl-47.4%  H2-50%  Ar  mixture  at  a post  shock  temperature  of  1070“K. 

Again,  the  solid  lines  are  computed  temporal  histories  for  various  assumed 
values  of  the  V-V  rate,  and  the  dashed  line  represents  the  relaxation  profile 
that  would  be  expected  for  a negligible  contribution  from  the  V-V  process. 

It  would  appear  that  the  experimental  data  obtained  from  this  run  would  best 
be  fit  with  a value  of  the  V-V  rate  constant  approximately  1/3  of  the  base- 
line value.  Analysis  of  other  shock  tube  data  taken  for  mixtures  of  high 
H2/HCI  concentrations  at  various  post  shock  temperatures  indicated  similar 
results;  that  is,  the  experimental  data  are  best  fit  with  a value  of  the  V-V 
rate  constant  between  1/3  to  3 times  the  theoretically  derived  baseline 
value.  Determination  of  the  V-V  rate  constant  to  better  precision  will  require; 
(a)  a more  precise  determination  of  the  HCI-H2  V-T  rate  constant  and,  (b) 
acquisition  of  more  highly  resolved  data  taken  for  higher  H2/HCI  concentration 
ratios . 

An  additional  test  may  be  applied  as  a constraint  on  the  value  of 
the  V-V  rate  constant.  Using  the  value  for  the  HCI-H2  V-T  rate  obtained  from 
the  present  experiments,  one  can  analyze  the  data  of  Bott  and  Cohen  near  800°k 
to  see  what  value  of  the  V-V  rate  constant  is  required  to  best  fit  their  data. 
Presented  in  Figure  2.7  are  the  results  of  such  an  analysis.  The  three  solid 
curves  shown  are  predicted  temporal  profiles  of  the  HCl  laser  Induced  fluores- 
cence decay  calculated  for  assumed  values  of  the  V-V  rate  constant  of  1/2,  1 
and  2 times  the  baseline  theoretical  value  at  770 °K.  The  dashed  curve  is  the 
relaxation  profile  observed  by  Bott.  It  is  apparent  that  the  data  can  best 
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Figure ,2.8  Semilogarithmic  Plot  of  Observed 
HCl  Vibrational  Relaxation  for 
Mixture  of  High  H2/HCI  Concen- 
tration Ratio  and  Comparison  with 
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be  fit  with  a value  for  the  V-V  rate  constant  very  near  the  baseline  value. 

A similar  comparison  made  with  Bott*s  data  taken  at  T ■ 554®K  Is  shown  In 

Figure  2.10  ^d  Indicates  essentially  the  same  result  - the  data  are  best 

fit  with  a value  for  the  V-V  rate  constant  very  near  our  previously  assumed 

theoretical  value.  In  addltloh,  a recent  as  yet  unpublished,  laser 

fluorescent  experiment  In  a HF-ici-H-  mixture  at  room  temperature  was  re- 
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ported  by  Bott  . Analysis  of  these  data  for  the  rate  of  the  H2/HCI  V-»V 
process  yields  a value  which  Is  within  70%  of  the  base  line  value  at  300®K. 
Therefore,  It  Is  concluded  that  the  present  experimental  data  support  our 
theoretically  derived  H2/HCI  V-+V  rate  constant  from  300  to  'v>1000°K,' 

Table  2.2  lists  the  current  values  of  rate  constants  for  the  vibrational 
energy  exchange  processes  for  the  H2/HCI  system. 

2.4  H-Atom  Kinetics 

In  the  Phase  I Investigation  the  potentially  serious  performance 
degradation  of  atoms,  particularly  the  H-atom,  produced  In  the  high  tempera- 
ture stagnation  region  and  subsequently  frozen  In  the  rapid  expansion,  was 

discussed.  At  that  time  only  one  measurement  of  HCl  (v  “ 1)  deactivation  by 

18 

H-atoms  was  available  , and  that  measurement  was  at  room  temperature.  Hence, 
a large  and  uncertain  temperature  extrapolation  was  necessary.  Table  2.3 
summarizes  the  previously  derived  rate  expressions  for  the  atomic  rate  pro- 
cesses. 

Early  In  the  current  program  an  unpublished  result  was  brought  to 

20 

our  attention*  that  the  rate  In  reference  18  had  been  confirmed  . Working 
over  the  limited  temperature  range  of  200  - 400 “K,  an  activation  energy  of 
about  1 kcal/mole  had  also  been  determined.  This  Information  has  been  con- 
verted Into  a temperature  dependent  rate  constant  and  Is  shown  on  Figure  2.11 
and  conq>ared  to  the  rate  based  on  Equation  (9)  of  Table  2.3.  Clearly,  the 
data  of  Glass  et  al,  lead  to  a steeper  temperature  dependence  and  larger 

rate  for  temperatures  above  lOO^K.  Also  shown  In  Figure  2.H  Is  the  theore- 

21 

tlcal  calculation  by  R.  Wllkens  , which  has  been  divided  by  10.  More 

19 

recently,  another  measurement  at  room  temperature  has  provided  a further 
confirmation  of  the  low  temperature  rate.  This  result  Is  also  shown  In  Figure  2.11. 

*Prlvate  communication  with  Dr.  Thomas  Barr  (USAMICOM) 
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Figure  2.10  Interpretation  of  Bott's  Observed 
HCI/H2  Data  in  Terms  of  V-V 
and  V-T  Processes 


TABLE  2.2 


SUMMARY  OF  UPDATED  VIBRATIONAL  RATE  CONSTANTS 
FOR  HARMONIC  OSCILLATOR  MODEL 


I 

N3 

I 


V ^ T Processes 

(1)  HCl  (v  = 1)  + HCl  (v  = 0)  ^ HCl  (v  = 0)  + HCl  (v  = 0) 

(2)  HCl  (v  = 1)  + H2  (v  = 0)  HCl  (v  - 0)  + H^  (v  = 0) 

(3)  H2  (v  = 1)  + H2  (v  = 0)  ^ H2  (v  = 0)  + H2  (v  - 0) 

(4)  H2  (v  = 1)  + HCl  (v  = 0)  ^ H2  (v  = 0)  + HCl  (v  “ 0) 


k (cm^  - part  ^ - sec 

- 4.68  X 10~^°  exp  (-76.8T~^^^) 

+ 3.13  X 10~^^  exp  (27.9 
k2  * 2.5  X 10"^^  exp  (-55.9 
k^  = 2.18  X 10"^  exp  (-114.9 
k^  - 4.35  X 10"^°  exp  (-114. 9T“^^^) 


V ^ V Process 


k = 1.71'^^;^^  X 10“^^  T exp  (-30.0 
^ -0.86 


(5) 


H2  9v  = 1)  + HCl  (v  = 0)  ^ H2(v  = 0)  + HCl  (v  = 1) 


TABLE  2.3 


VIBRATION  RATE  CONSTANTS  FOR  ATOMIC  SPECIES  FOR 

H^/HCl  SYSTEM 

3 -1  -1 

(Rate  constants  in  units  of  cm  - particle  - sec  ) 

(Rate  expressions  are  for  the  forward  process) 


V ->■  T Processes 


HCl  (v  = 1)  + H ^ 
k = 3.75  X 10 


HCl  (v 
-13  ^1/2 


0)  + H 


(9) 


HCl  (v  = 1)  + Cl  ^ HCl  (v  = 0)  + Cl 
k = 7.50  X 10~^^  T^'^^ 


(17) 


H^  (v  = 1)  + H ^ H2  (v  - 0)  + H 

I 1 in  ia“1^  1/2 

k = 1.73  X 10  T 


(18) 


H2  (v  = 1)  + Cl  2 ^2 

-15  1/2 

k = 1.73  X 10  t' 


(19) 
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Figure  2.  11  Rate  of  Deactivation  of  HCl  (v  = 1 
by  H-atoms.  Comparison  of 
Experimental  Data,  Temperature 
Extrapolations,  and  Theory 


Clearly,  the  data  at  and  around  roo.^  teinpp’"ature  are  in  close  agree- 
ment and  support  a large  rate  for  the  deactivation  of  HCl  (v  = 1)  by  H-atoms. 

The  discrepancy  between  the  measured  values  for  HCl/H  and  HF/H  still  exists 
and  requires  adequate  theoretical  explanation.  Nevertheless,  the  large  room 
temperature  value  of  HCl/H  combined  with  the  uncertainty  In  temperature  depen- 
dence make  this  process  of  key  Importance  in  understanding  potential  laser 
performance.  No  additional  data  on  the  rates  of  other  processes  shown  In 
Table  2.2  or  2.3,  have  become  available  during  this  program. 

2.5  H-HCl  Deactivation  Experiment  Design 

An  experiment  to  measure  the  rate  of  H-atom  deactivation  of  HCl  for 
T > 1000°K  la  shown  schematically  in  Figure  2.12.  The  basic  concept  consists 
In  converting  the  present  1.5  Inch  diameter  shock  tube  into  a glass  flow 
reactor/shock  tube.  A primary  flow  of  H2  and  Ar  will  be  passed  through  a micro- 
wave  discharge  and  mixed  upstream  with  a secondary  flow  of  HCl  external  to  the 
shock  tube.  This  mixture  will  then  pass  through  the  shock  tube.  Flow  rates 

of  H , Ar,  and  HCl  will  be  measured  ttslng  mechanical  flow  meters.  H-atom 
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concentrations  will  be  determined  by  the  standard  NOCl  titration  technique 

Once  a steady  flow  of  determined  concentration  is  established,  fast  acting 

valves  will  close  off  the  tube,  and  a shock  wave  will  be  initiated.  I.R. 

emission  from  the  HCl  will  then  be  monitored  behind  the  Incident  shock  wave 

using  the  InSb  detector  as  In  the  present  experiment 

In  designing  the  above  experiment,  various  prereactions  that  might 

occur  in  the  flow  must  be  considered.  These  chemical  reactions  are: 


H + HCl  H2  + Cl 

(10) 

Cl  1/2  Clj 

(11) 

H 1/2 

(12) 

H + CI2  ^ HCl  + Cl 

(13) 

H + H+  M-J-H2+M 

(U) 

Cl  + Cl  + M CI2  + M 

(15) 

H + Cl  + M HCl  + M 

(16) 
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Figure  2.12  Shock  Tube  Flow  Reactor  for 

H-Atom  Deactivation  Experiment 


Rate  constants  for  these  various  reactions  were  obtained  and  a computer  pro- 
gram used  to  calculate  the  flow  history  of  the  relevant  species.  Table  2.4 
lists  the  values  of  the  rate  constants  utilized  in  the  program.  For  the 
flow  pressures  of  lnterest2^^2^all  accomodation  coefficients,  y>  for  H and 
Cl  atoms  of  less  than  10  ^ , the  rate  limiting  processes  for  wall. removal 

of  atoms  are  calculated  to  be  reactions  (11)  and  (12)  and  not  diffusion.  See 
reference  (24)  for  a fuller  explanation  of  the  theory  of  wall  effects  on  active 
species . 

One  of  the  principal  results  of  the  above  flow  kinetics  calculations 
was  that  even  at  the  lowest  Aci  concentrations  permissible  for  adequate 
fluorescence  signal  and  for  the  fastest  flow  speeds  reasonably  attainable 
through  the  shock  tube,  the  time  constant  associated  with  reaction  (10)  Is 
expected  to  be  short  compared  to  the  flow  time  between  the  HCl  injector  and 
the  teat  region.  Such  a situation  Is  undesirable  because  It  produces  both 
a depletion  of  H-atoms  and  an  undesirable  concentration  of  Cl  atoms.  Con- 
cluding that  It  la  not  possible  to  design  a practical  shock  tube  flow  exper- 
iment to  compete  with  the  rapid  rate  of  reaction  (10) , it  was  decided  to 
prevent  reaction  (10)  from  proceeding  too  far  in  the  forward  direction  by 
operating  under  conditions  where  the  H2  concentration  Is  large  compared  to 
HCl,  and  the  backward  reaction  Is  favored.  Shown  in  Figure  2.13  are  the 

results  of  a calculation  for  a pre-shock  pressure  of  10  torr,  an  ambient 

3 

temperature  of  300®K,  a flow  velocity  of  10  cm/sec,  and  a mixture  of  1%  HCl  - 

33%  H2  - 0.3%  H - 65.7%  Ar.  Plotted  are  the  concentration  profiles  as  a 

function  of  flow  distance  down  the  tube.  As  can  be  seen,  for  a wall  accomo- 

-5  -4 

dation  coefficient,  y = 10  or  10  (values  that  can  easily  be  achieved  on 
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a properly  treated  glass  wall)  , the  H-atom  concentration  will  remain 
essentially  constant,  and  there  will  be  negligible  formation  of  Cl-atoms. 

Operating  In  mixtures  of  high  H2/HCI  ratios  does,  however,  present 
an  additional  potential  problem.  The  energy  sharing  between  H2  and  HCl  due 
to  the  V-V  process  can  complicate  the  Interpretation  of  the  effect  of  H-atom 
deactivation.  That  is,  if  HCl  and  H2  are  strongly  coupled  together,  does 
one  have  to  consider  the  contribution  from  H-atom  deactivation  of  H2 , as 
well?  Figure  2.14  shows  the  results  of  a calculation  which  Indicates  an 
experimental  regime  where  the  Interpretation  of  the  observed  HCl  relaxation 
is  essentially  Independent  of  the  effect  of  the  H on  H2.  The  calculation 
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TABLE  2.4 


LIST  OF  RELEVANT  RATE  CONSTANTS  FOR  H-ATOM  FLOW  EXPERIMENT 


H + HCl  ^ H2  + Cl  (10) 

k^o  " 5.8  + 2.5  X lO"^^  T^^^  exp  ((-1450  + 150) /T] 

Cl  1/2  CI2 

H 1/2  H2 

H + Cl„  ^ HCl  + Cl  (13) 

z ^ 

" 7 X 10"^°  exp  (-1500/T) 


k ■ i^ere  y is  the  wall  accomodation 
coefrlclent,  "c  la  the  mean  velocity  of  recom- 
hlnlng  species,  and  r Is  the  tube  radius. 


(11) 


(12) 


H + 


H + M H2  + M 
k,,  - 4.5  X 10 


(14) 


Cl  + Cl  + M CI2  + M (15) 

-29  -1 

k^^  - 3 X 10  T 

i 

H + Cl  + M HCl  + M (16) 

k^g  - 4 X 10"^®  t“^ 
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Ref. 

22.23 

22 

22.24 

25 

26 
25 
25 


SPECIES  MOLE  FRACTION 


P*  lOtorr 
T«300 

V = l0^  cm/sec 
[HCI]  = .01 
« .33 

CAr3  = .657 


0 


10 


Figure  2.13  Results  of  Flow  Tube  Kinetics 

Calculation  for  Design  of  H-Atom 

Deactivation  Experiment 


0 1.0  2.0  3.0  4.0  5.0  6.0 

^lob> 

Figure  24.4  Shock  Tube  Relaxation  Calculation 
with  H-Atoms  Showing  Decoupling 

of  HCl  and  H2  Relaxation 
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was  performed  for  an  incident  shock  speed  of  1.35  mm/ysec  In  a 5 torr  mixture 
of  2.5%  HCl  - 4.5%  H - 43%  H2  - 50%  Ar.  Such  conditions  yield  a post  shock 
temperature  of  1030°K.  Even  at  this  Increased  H-atom  concentration,  calcula- 
tions indicate  that  3 body  recombination  of  H is  still  slow  compared  to  the 
pertinent  flow  time.  The  two  dash-dot  curves  plotted  in  Figure  2.13  show  the 

predicted  relaxation  profile  of  the  H.  vibrational  energy  assuming  our  present 

1 ^ 

base  line  value  for  the  H2/H  rate  (upper  curve) , and  a rate  ten  times  that 
value  (lower  curve) . The  two  solid  curves  plotted  show  the  predicted  relaxa- 
tion of  the  HCl  vibrational  energy  for,  again,  an  assumed  value  of  the  H2/H 
deactivation  rate  equal  to  our  baseline  rate,  and  10  times  that  value. 

Clearly,  the  value  for  the  rate  of  H-atom  deactivation  of  H2  has  little 
effect  on  the  HCl  relaxation  profile  under  these  conditions.  Similar  calcu- 
lations were  performed  for  other  conditions,  and  the  conclxision  is  reached 
that  a determination  of  the  rate  for  H-atom  deactivation  of  HCl  can  be  made 
independent  of  the  rate  for  H-atom  deactivation  of  H2  provided  that  the 
energy  sharing  between  H2  and  HCl  is  not  allowed  to  progress  appreciably. 

Such  a condition  can  be  achieved  by  operating  at  the  lower  shock  temperatures 
and/or  creating  sufficient  H-atom  loading  to  make  the  V-T  processes  fast 
compared  to  the  rate  for  V-V  coupling. 

2.6  HCl  Spectral  Characteristics 

Besides  the  basic  kinetic  processes  discussed  above,  it  was  decided 

to  assess  two  other  important  parameters  used  in  the  modeling  - the  Einstein 

coefficients  and  line  broadening  cross-sections  for  HCl.  Both  of  these  latter 

parameters  are  needed  in  the  calculation  of  gain  and  subsequent  power  extraction. 

The  experimental  data  on  Einstein  coefficients  for  the  fundamental 

vibration/rotation  band  of  HCl  exist  for  only  the  lower  vibrational  levels. 

Theory  must  be  used  to  extrapolate  these  data  to  the  upper  levels.  The  results 

28  29 

of  two  such  theoretical  extrapolations  ’ are  shown  in  Figure  2.15.  These 

results  are  in  close  agreement  at  the  lower  vibrational  levels  but  exhibit  an 

increasing  divergence  with  increasing  v.  While  this  descrepancy  is  of  some 

fundamental  importance,  its  impact  on  the  present  laser  modeling  is  small, 

since  the  emphasis  is  on  energy  extraction  from  the  lower  vibrational  levels 

to  provide  laser  wavelengths  with  good  atmospheric  transmission.  For  the  pre- 

28 

sent  modeling,  the  Einstein  coefficients  of  Smith  are  used. 
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Figure  2.15  Einstein  Coefficient  for  HCl 
v-»v-l  transitions 
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Another  important  physical  parameter  required  in  the  laser  modeling 
is  the  optical  broadening  cross-section.  The  experimental  data  base  for  HCl 

is  somewhat  limited  and  is  shown  in  Figure  2.16  as  the  half  width,  y,  for  the 

30  31  , 

P-branch  v " 1 -+■  0 transition  * . The  collision  broadening  shows  an  interest- 

ing variation  with  rotational  level.  This  type  of  variation  has  also  been 

32 

noted  with  other  hydrogen  halides  such  as  HF  . Also,  the  self-broadening 
(HCl  by  HCl)  is  considerably  larger  than  other  homonuclear  species  such  as  H2 . 
Finally,  the  data  of  these  two  investigations  for  similar  species  are  in  Close 
agreement.  No  data  exist  for  the  higher  vibrational  levels  of  HCl  and  for 
higher  temperature. 

Meredith  and  co-workers  have  extensively  modeled  the  optical  broad- 
ening cross-sections  for  HF  and  DF  in  response  to  the  needs  of  chemical,  laser 

32 

development.  For  example,  they  have  shown  that  the  cross-section  can  vary 
as  T to  T depending  upon  the  type  of  molecular  interaction  involved. 

Lacking  this  type  of  theoretical  calculation  for  HCl,  we  have  determined  an 
optical  broadening  cross-section  from  the  data  base  of  Figure  2.16  and 
assumed  no  temperature  or  vibrational  level  dependence.  For  high  J transitions 
(J  > 8) , we  have  assumed  that  the  cross-section  is  constant  and  equal  to  its 
value  at  J = 8.  Clearly,  refinement  of  those  assumptions  is  necessary  as  the 
concept  of  an  H2/HCI  GDL  advances. 
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3.0 


PROPELLANT  DEVELOPMENT 


3.1  Objective 

The  objective  of  the  propellant  development  is  to  develop  solid  propel 
lant  gas  generators  that  yield  combustion  products  suitable  for  use  as.  the 
working  fluid  in  an  advanced  H2-HCI  gas  dynamic  laser.  There  are  certain  basic 
constraints  concerning  the  properties  of  the  combustion  products  produced  by 
the  gas  generators.  These  constraints  are  enumerated  below: 

1.  Flame  Temperature 

Initially,  in  Phase  I it  was  believed  that  the  required  flame 
temperature  value  was  3000°K,  but  nozzle  optimization  calculations  for  available 
power  indicated  that  a decrease  in  flame  temperature  was  permlssable.  In  addi- 
tion, at  lower  temperatures  there  are  less  H-atoms  present.  As  shown  ih  Figure 
3.1,  H-atoms  increase  with  temperature  due  to  H2  dissociation.  Thus,  at  a 
lower  temperature  the  H-atom  problem  is  somewhat  alleviated.  A parametric  study 
of  available  power  for  a range  of  flame  temperatures  in  the  range  2400  to  3000 °K 
indicated  2800°K  to  be  the  beat  compromise  value. 

2 . Combustion  Pressure 

A chamber  or  combustion  pressure  of  750  psia  was  selected.  This 
pressure,  and  the  above  flame  temperature,  are  the  stagnation  values  for  the 
working  fluid. 

3.  Gas  Composition 

The  molar  ratio  selected  (H2/HCI)  in  the  combustion  products  is 
97/3.  Initially,  values  as  low  as  90/10  were  considered.  More  recent  results 
indicate  that  an  optimum  (H^/HCl)  value  is  97/3  and  values  as  hi^  as  98/2  may 
be  desirable. 

4.  Gas  Purity  (Molar  Content) 

The  H2/HCI  purity  in  the  combustion  products  must  be  >90%. 
Initially,  the  requirement  was  for  a best  possible  value,  and  hydrocarbons  - 
which  are  an  excellent  off-the-shelf  source  of  hydrogen  - were  considered 
seriously  as  propellant  ingredients.  However,  subsequent  calculations  of  the 
effect  of  contaminants  primarily  CO  and  N2  on  gain  coefficient  and  available 
power  indicated  that  a relatively  pure  gas  composition  for  the  combustion 
products  is  required. 
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5.  Condensed  Species 

The  only  way  to  generate  the  flame  temperatures  of  Interest  Is 
by  the  use  of  tnetal  fuels  (principally  A1  and  Zr) • Hydrocarbons  are  not  good 
fuels  because  combustion  to  CO  Is  not  energetic  enough  while  combustion  In 
excess  of  CO  yields  H2O,  which  is  detrimental  to  laser  performance.  Combustion 
of  A1  and  Zr  results  In  the  production  of  AI2O2  and  Zr02»  both  of  which  are 
condensed  at  the  flame  temperatures  and  pressure  of  Interest,  and  must  be  re- 
moved from  the  combustion  products,  either  by  clinker  formation  or  by  mechanical 
means  (filtration  and/or  inertial  separation). 

6.  Weight  Yield 

The  weight  yield  of  the  H2/HCI  mixture  should  be  as  high  as 
possible  by  weight  of  propellant.  The  results  of  a preliminary  laser  systems 
analysis  Indicates  that  a 5%  yield  Is  desirable.  This  5%  yield,  however,  was 
not  considered  to  be  a hard  requirement,  and  two  propellants  were  developed 
which  were  satisfactory  and  feasible  In  all  respects  except  with  respect  to 
weight  yield. 

The  one  propellant  ingredient  that  Is  necessary  to  meet  all  of  the 
above  requirements  Is  hydrazine  dlborane  referred  to  as  HDB  through- 

out this  report) . The  compound  Is  available  in  limited  quantities  on  a 
special  order  basis,  but  synthesis  is  a straightforward  and  a well-known  process. 

Three  candidate  propellant  compositions  were  developed  under  this 
program.  The  three  are  pressed  compositions  comprised  of  an  oxidizer,  a fuel, 
a coolant  and  a chloride.  Candidates  are  based  upon  selection  of  ingredients  and 
computed  values  of  theoretical  flame  temperature  and  combustion  product  composi- 
tion. For  a given  set  of  ingredients,  propellant  composition  Is  systematically 
varied  to  map  (H2  + HCl)  weight  yield  and  purity  versus  the  molar  ratio  (H2/HCI) 
within  the  constraint  T^  =•  2800'’K.  Candidate  compositions  were  then  charac- 
terized for  safety  and  ballistic  properties  and  fired  in  subscale  motors. 
Additionally,  various  filter  arrangements  were  evaluated  to  retain  condensed 
combustion  products  in  the  motor.  The  methodology  and  results  of  this  work  are 
described  below. 


- 36  - 


, , TABLE  3.1 

EQUILIBRIUM  YIELD  VALUES  AND  HEATING  VALUES  0 P =750  PSIA  FOR  SELECTED  PROPELLANT  INGREDIENTS 


Equilibritim  Yields, 
mol/100  ems  0 2800“K 

PROPELLANT  COMPOSITION,  WT  % 

«2 

HCl 

2ng 

40.467%  N2H^N02  + 59.533%  ZrH2 

36.402%  NH,N0„  + 63.598%  ZrH- 
4 3 2 

38.654%  NH.CIO,  + 61.346%  ZrH„ 

4 4 2 

’63.078%  Zr(OH)^  + 36.922%  ZrH2 
45.025%  Zr0(N02)2*2H20  + 54.975%  ZrH2 


NH,N0„  + 63.598%  ZrH_ 

4 3 2 

NH.CIO,  + 61.346%  ZrH„ 
4 4 2 

Zr(OH)^  + 36.922%  ZrH2 


38.041%  Zr(N02)^*5H20  + 61.959%  ZrH2 

62.375%  WO^  + 37.625%  ZrH2 

57.264%  H.WO,  + 42.726%  ZrH. 

2 4 2 

59.972%  NC  + 40.028%  ZrH 


57.776%  TMETN  + 42.224%  ZrH, 


69.064%  NitroGuanidine  + 30.936%  ZrH2 
32.783%  A1(N(^)2'9H20  + 67.217%  ZrH2 
52.730%  Al(OH).,  + 47.270%  ZrH., 


AKOH)^  + 47.270%  ZrH2 


* 'See  Table  3.2  for  Glossary  of  Terms 


1.685 

1.576 


1.181 

0.9180 

1.098 

0.4012 

0.6835 

1.241 

1.450 

1.635 

1.492 

1.506 


2.351 

2.055 


1.144  0.3048  1.649 


Equilibrium 
K cal/100 


1.195 
1.100 
1.291  * 
0.4059 
0.6^15 
2.855 
2.948 
3.657 
1.650 
1.533 


3000  2800  2600 


-19.98  -40.91  -47.98 

-18.77  -40.13  -46.78 

-47.50  -67.44  -73.25 

52.70  30.48  25.17 

-19.74  -40.82  -45.78 

-33.73  -55.14  -60.48 

6.762  - 4.599  - 7.347 

11.74  - 1.565.  - 5.030 

1.287  -15.40  -22.58 

- 8.813  -26.35  -33.83 

28.75  12.68  4.269 

- 6.505  -29.01 

67.03  47.67 


Heating  Values, 
Sj  0 T.  °K 


2600  2400 


-54.66 

-53.06 

-78.74 

20.11 


3.2 


Selection  of  Propellant  Ingredients 

The  propellants  are  comprised  of  oxidizer,  fuel,  coolant  and  chloride, 
each  Ingredient  selected  on  the  basis  that  the  constraints  of  the  problem  are 
achieved.  Candidate  Ingredients  were  selected  on  the  basis  of  hydrogen  content, 
chemical  energy  content,  commercial  availability,  ease  of  use,  stability  of 
propellant  Ingredients  and  stability  of  combustion  products. 

For  any  given  set  of  propellant  Ingredients,  there  is  a unique  compo- 
sition that  yields  the  desired  values  of  flame  temperature  and  H2/HCI  composi- 
tion. To  minimize  the  appearance  of  H^O  as  a combustion  product,  the  stoichio- 
metry of  candidate  propellants  was  constrained  to  an  oxidation  ratio  (0  ) of 

R 

unity,  l.e.,  the  oxygen  in  the  propellant  Is  equal  to  that  required  to  burn  the 
metals  to  the  metal  oxides  (Zr02,Al202)  and  the  carbon  (if  present)  to  CO.  The 
oxides  must  be  stable  with  respect  to  at  the  stagnation  or  chamber  conditions 
to  prevent  formation  of  H2O.  The  resultant  stoichiometry  is  referred  to  through- 
out this  report  as  the  OMOX  point.  Compositions  fuel- rich  of  this  stoichiometry 
are  not  particularly  desirable  because  metal  vapors  and  suboxides  are  produced 
as  combustion  products,  which  may  interfere  with  lasing  or  result  In  condensed 
species  in  the  optical  cavity. 

Equilibrium  values  of  H2  yield,  HCl  yield,  total  gas  yield  (Eng), 
and  heating  value  (AH°)  were  computed  for  selected  OMOX  combustions  of: 

(1)  oxidizer  and  fuel;  (2)  coolant  and  either  oxidizer  or  fuel,  and  (3)  chloride 
(and  fuel  if  necessary) . Calculations  were  carried  out  with  the  ARC  equili- 
brium program  using  JANAF  thermo chemical  data  on  an  IBM  370/168  computer  at 
assigned  values  of  750  psia  combustion  pressure  and  temperatures  of  3000  to  2400°K 
with  200 °K  increments.  The  results  of  these  calculations  are  given  in  Table  3.1. 
The  heating  values  in  the  table  are  derived  quantities  which  are  defined  as  the 
enthalpy  at  equilibrium  of  the  combustion  products  at  a temperature  T and  750  psia 
less  the  enthalpy  of  the  propellant.  Negative  values  of  AH°  indicate  available 
heat  while  positive  values  indicate  that  heat  is  required  to  attain  the  tempera- 
ture T.  For  example,  the  OMOX  composition 

N2H^N02  + 1.5  ZrH2  1.5  Zr02  + 4H2  + 1.5  N2 
is  comprised  of  40.467%  N2H^N02  + 59.533%  ZrH2,  (by  weight)  and  has  heating 
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TABLE  3.1 
(Continued^ 

EQUILIBRIUM  YIELD  VALUES  AND  HEATING  VALUES  @ P =750  PSIA  FOR  SELECTED  PROPELLANT  INGREDIENTS 

c 


PROPELLANT  COMPOSITION,  WT  % 

Equilibrium  Yields, 
mol/100  ms  @ 2800“K 

Equilibrium  Heating  Values , 
K cal/100  gms,  @ T,  °K 

«2 

HCl 

Ing 

3000 

2800 

2600 

2400 

40.995%  N2H^N02  + 59.005%  Zr 

1.066 

- 

1.732 

-59.96 

-79.64 

-85.44 

-90.99 

36.909%  NH.NO,  + 63.091%  Zr 
4 3 

0.9124 

- 

1.389 

-61.49 

-86.81 

-86.81 

-91.88 

39.174%  NH.CIO,  + 60.826%  Zr 
4 . 4 

0.4972 

0.3051 

1.000 

-89.06 

-107.7 

-112.2 

-116.5 

63.581%  Zr(OH)^  + 36.419%  Zr 

0.7938 

- 

0 . 8031 

28.625 

7.156 

- 1.699 

62.886%  WO^  + 37.114%  Zr 

0.0 

- 

0.0 

-18.15 

-28.64 

-30.56 

-32.45 

57.798%  H-WO,  + 42.202%  Zr 
2 4 

0.2300 

- 

0.2326 

-16.54 

-28.87 

-31.40 

-33.87 

60.495%  NC  + 39.505%  Zr 

0.8233 

- 

2.446 

-25.32 

-41.10 

-47.42 

-53.49 

58.308%  TMETN  + 41.692%  Zr 

1.011 

- 

2.517 

-36.96 

-53.56 

-60.15 

-66.45 

69.529%  Nitro^uanldlne  + 30.471%  Zr 

1.315 

- 

3.346 

8.407 

- 6.956 

-14.72 

-22.09 

63.791%  N2H^N02  + 36.209%  A1 

1.657 

- 

2.700 

-89.38 

-101.6 

-112.8 

-123.4 

59.734%  NH.NO-  + 40.266%  Zl 
4 3. 

1.475 

• - 

2.252 

-95.52 

-107.2 

-118.1 

-128.5 

62.021%  NH.CIO,  + 37.979%  A1 
4 4 

0.7689 

- 

1.593 

-135.2 

-145 . 7 

-155.5 

-164.7 

74.301%  AKOH)^  + 25.699%  A1 

1.415 

- 

1.440 

51.79 

40.20 

29.30 

18.84 

77.643%  H-WO,  + 22.357%  A1 
2 4 

0.3077 

- 

0.3131 

-20.47 

-25.24 

.-29.82 

-34.26 

81.573%  Zr(OH)^  + 18.427%  A1 

1.014 

1.032 

38.76 

20.63 

.13.72 

7.126 

TABLE  3.1 
(Continued> 

EQUILIBRIUM  YIELD  VALUES  AND  HEATING  VALUES  @ = 750  PSIA  FOR  SELECTED  PROPELLANT  INGREDIENTS 


PROPELLANT  COMPOSITION,  WT  % 

Equilibriiim  Yields, 
mol/100  ems  @ 2800®K 

Equilibrium  Heating  Values, 
K cal/100  gms,  @ T,  °K 

«2 

HCl 

Ing 

3000 

2800 

2600 

2400 

69.316%  N2H^N02  + 30.684%  CH^ 

3.955 

- 

7.305 

163.8 

147.8 

133.2 

119.5 

65.544%  NH.NO,  + 34.456%  CH. 

4 3 2 

4.038 

- 

7.383 

188.7 

172.4 

157.7 

143.8 

67.680%  NH.CIO,  + 32.320%  CH. 

4 4 2 

3.129 

0.5697 

6.353 

125.5 

111.4 

98.54 

86.61 

85.022%  Zr(OH)^  + 14.978%  CH2 

2.114 

- 

3.214 

167.9 

147.5 

139.1 

131.1 

84.638%  WO^  + 15.362%  CH2 

1.082 

- 

2.196 

108.4 

102.9 

97.92 

93.21 

81.662%  H-WO,  + 18.338%  CH» 
2 4 2 

1.616 

- 

2.950 

135.2 

128.1 

121.6 

115.5 

73.134%  Zr0(N02)2*2H20  +26.866%  CH2 

2.428 

- 

4.660 

141.5 

124.6 

114.3 

67.112%  Zr(N02)^*5H20  + 32.888%  CH2 

3.074 

5.804 

154.0 

136.4 

124.2 

68.857%  N2H^N03  + 31.143%  HTPB 

3.414 

6.796 

138.0 

123.1 

109.6 

96.92 

65.056%  NH.NO.  + 34.944%  HTPB 
4 3 • 

3.430 

- 

6.812 

159.4 

144.6 

131.0 

118.3 

67.208%  NH.CIO,  + 32.792%  HTPB 
4 ‘4 

2.565 

0.5654 

5.824 

98.54. 

85.75 

' 

73.96 

63.03 

TABLE  3.1 
(Continued) 

EQUILIBRIUM  YIELD  VALUES  AND  HEATING  VALUES  @ P = 75Q  PSIA  FOR  SELECTED  PROPELLANT  INGREDIENTS 


PROPELLANT  COMPOSITION,  WT  % 

Equilibrium  Yields, 
mol/100  gms  @ 2800“K 

Equilibrium  Heating  Values, 
K cal/100  gms,  @ T,  “K 

«2 

HCl 

£ng 

3000 

2800 

2600 

2400 

29.465%  NH.NO-  + 70.525%  PEG  200 
4 3 

3.949 

7.299 

183.5 

167.4 

152.8 

139.1 

31.49^3%  NH.CIO,  + 68.507%  PEG  200 
4 4 

3.530 

• 

0.2652 

6.823 

154.2 

139.2 

125.4 

112.6 

52.223%  NC  + 47.777%  PEG  200 

2.894 

- 

6.306 

135.7 

122.0 

109.5 

97.77 

'49.957%  TMETN  + 50.043%  PEG  200 

3.153 

- 

6.550 

133.3 

119.1 

106.1 

93.84 

61.959%  Nitroguanidine  + 38.041%  PEG  200 

2.907 

- 

6.319 

132.9 

119.2 

106.7 

94.95 

55.552%  NH.NO-  + 44.448%  TCE 
4 3 

1.360 

— 

4.863 

63.21 

53.10 

43.68 

34.73 

57.910%  NH.CIO,  + 42.090%  TCE 

0.4838 

4.115 

45.22 

36.49 

28.41 

20.77 

Cf 

76.589%  NC  + 23.411%  TCE 

1.038 

- 

4.558 

49.93 

40.58. 

31.81 

23.44 

74.923%  TMETN  + 25.077%  TCE 

1.294 

- 

4.800 

40.37 

30.40 

21.11 

12.2  7 

100%  NH.Cl 
4 

2.791 

1.849 

5.641 

227.6 

214.4 

202.7 

192.0 

100%  N2H^C1 

2.903 

1.444 

5.870 

169.1 

155.6 

100%  HDB 

8.200 

■ 

8.546 

227.7 

128.9 

87.65 

59.99 

EQUILIBRIUM  YIELD  VALUES  AND 


PROPELLANT  COMPOSITION,  WT  % 

79.923%  Zr(OH)^  + 20.077%  AlH^ 

57.154%  NH^NO^  + 42.846%  AlH^ 

79.439%  WO^  + 20.561%  AlH^ 

75.745%  H-WO,  + 24.255%  AlH, 

2 4 3 

61.303%  N2H^N0^  + 38.697%  AlH^ 
59.489%  NH^CIO^  + 40.511%  AlH^ 
72.221%  A1(0H)2  + 27.779%  AlH^ 


TABLE  3.1 
(Continued) 

; VALUES  @ P =750  PSIA  FOR  SELECTED  PROPELLANT  INGREDIENTS 
c 


Equilibrium  Yields, 
mol/lOO  gms  @ 2800°K 


Equilibrium  Heating  Values, 

1 K cal/100  gms,  @ T, 

“K 

3000 

2800 

2600 

2400 

HCl 


%ng 


1.988 

3.532 

1.018 

1.501 

3.508 

2.727 

2.750 


0.4244 


2.023 

4.313 

1.036 

1.527 

4.544 

3.578 

2.799 


64.16 

-37.62 

4.110 

10.45 

-37.31 

-77.42 

85.18 


42.60 
-54.82 
- 2.143 
2.422 
-54.44 
-93.34 
70.05 


34.52 

-70.16 

- 7.844 

- 4.823 
-69.68 
107.74 

56.26 


26.07 
-84.34 
-13.19 
-11.56 
-83.  75 
120.99 
43.3- 


TABLE  3.2 

GLOSSARY  OF  CANDIDATE  PROPELLANT  INGREDIENTS 

Chemical 
Formula  or 

INGREDIENT Acronym 

A.  OXIDIZER 


1. 

Hydrazine  Nitrate 

N2H5NO3 

2. 

Ammonium  Nitrate 

NH.NO. 
4 3 

3. 

Ammonium  Perchlorate 

NH  CIO 

4. 

Tungsten  Oxide 

5. 

Tungstic  Acid 

H2WO4 

6. 

Zirconyl  Nitrate  Dihydrate 

Zr0(N03)2*2H20 

7. 

Zirconium  (IV)  Nitrate  Pentahydrate 

Zr(N03)^-5H20 

8. 

Nitrocellulose  (A  Resin  used  in  Solid 
Propellant  Binders) 

NC 

9. 

Trlmethylole thane  Trinitrate  (A  co- 
plasticizer for  Nitrocellulose)' 

TMETN 

FUELS 

1. 

Zirconium 

Zr 

2. 

Zirconium  Hydride 

ZrH2 

3. 

Aluminum 

A1 

4. 

Aluminum  Hydride 

AIH3 

5. 

Tetracyanoethylene  (A  High-energy 
Organic  Compound,  C2(CN)^) 

TCE 

COOLANTS 

1. 

Zirconium  (IV)  Hydroxide 

Zr(OH)^ 

2. 

Aluminum  Hydroxide 

A1(0H)3 

3. 

Hydrazine  Diborane  (®2^2®10^ 

HDB 

4. 

Polyethylene,  Wax,  or  Long-chain 
Hydrocarbon 

CH2 
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TABLE  3.2 

(GLOSSARY  OF  CANDIDATE  PROPELLANT  INGREDIENTS  - CONTINUED) 


Chemical 
Formula  or 


5.  Hydroxy  Terminated  Polybutadiene  (An  HTPB 

Isocyanate-cured  Polyurethane  Binder 
System  Widely  used  In  the  Propellant 
Industry) 


6 . Polyethylene  Glycol  with  a Molecular 

Weight  of  200  a.m.u. 

PEG-200 

7 . N1 tro  guanidine 

CHLORIDES 

1.  Ammonium  Chloride 

NH^Cl 

2.  Hydrazine  Chloride 

N2H3CI 

3.  Zirconium  (IV)  Chloride 

ZrCl^’ 
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values,  Ah“  = -19.98,  -40.91,  -47.98  and  -54.66  JCcal/100  gm  @ 3000,  2800, 

2600  and  2400“K,  respectively.  These  heating  values  Indicate  the  heat 
available  at  the  Indicated  temperature.  Similarly,  the  OMOX  composition 

Zr(OH)^  + ZrH2  ^ 2 Zr02  + 3 H2 

Is  comprised  of  63.078%  Zr(OH)^  + 36.922%  ZrH2  (by  weight)  and  has  heating 
values  of  52.70,  30.48,  25.17  and  20.17  Kcal/100  gmat  3000,  2800,  2600  and 
2400°K,  respectively.  Indicating  that  heat  Is  required  to  attain  these  temper- 

O 

atures.  Note  that  in  the  above  two  examples  the  discontinuity  in  AH^Vs  temper- 
ature  between  2800®K  and  3000®K.  This  Is  due  to  the  heat  of  fusion  of  ZrO^ , which 
takes  place  at  2950“K.  Note  that  the  tabulated  values  of  H2,HC1  and  Eng  are  eaulll- 
-brlum  values  at  2800®K,  not  the  stoichiometric  values  expressed  above.  The  two 
do  not  differ  much,  especially  at  the  lower  temperatures  and  differ  primarily 
due  to  H2  dissociation. 

The  results  In  Table  3.1  are  discussed  In  terms  of  H2  yield,  H2 
yield  purity  (H2/Sng) , heat  release  and  oxygen  content  (of  oxidizers).  At 
the  start  of  the  program,  H2  yield  was  held  most  Important  and  the  effort 
was  directed  primarily  toward  maximizing  this  parameter.  For  this  reason 
organic  propellant  Ingredients  were  Initially  favored.  Later  on  In  the  program. 

It  became  apparent  that,  contaminants  such  as  CO  greatly  reduce  lasing  performance, 
and  for  this  reason  the  effort  was  redirected  toward  maximizing  H2  yield 
at  values  of  H2  yield  purity  >95%  and  the  use  of  organics  was  eliminated. 

Heat  release  and  oxygen  content  of  Oxidizers  are  most  Important  when  hlgh- 
H2~content  compounds  such  as  HDB  or  AlH^  are  used  In  the  propellants.  Salient 
features  of  the  results  In  Table  3.1  are  discussed  below. 

3.2.1  Oxidizers 

In  terms  of  H2  yield  the  Ingredients  that  generate  the  most  H2  In 

descending  order  are  N^H^NO.,  > KH^NO^  > NH^CIO^  > H2W0^  > WO^.  In  addition  to 

the  above  there  are  other  oxidizers  which  are  less  effective  than  NH,N0_  and 

4 3 

are  valued  In  the  order  AKNO^) 3*9820  > TMETN  > NC  > NH^ClO^  > Zr  (803)^ -5820  > 

ZrO (803)2 *2820  > 8280^.  In  terms  of  gas  purity  of  82  yield  (H2/Eng) , 8280^  = W03> 
NH^CIO^  > NH^N03  ^ ^2**5^°3’  heating  value,  NH^CIO^  > 88^803  > 

N2H5NO3  > WO3  > 8280^,  and  In  terms  of  oxygen  content,  88^803  > NH^CIO^  > 

N2H3NO3  > 8280^  > WO3. 
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All  of  these  oxidizers  are  commercially  available  and  relatively  easy 
to  use,  although  both  non-metalllc  nitrates  require  specially-dried  processing 
facilities  (relative  humidity  <io%  at  70**F)  for  use.  Of  the  two  nitrates,  N2H^N02 
Is  preferred  because  this  compound  does  not  undergo  a low-temperature  phase 
transition  as  does  NH^NO^.  Both  of  the  nitrates,  as  well  as  NH^CIO^  yield  sub- 
stantial amounts  of  N2  as  a combustion  product,  which  Is  a contaminant.  This 
problem,  however.  Is  not  severe  with  NH^CIO^  because  typically  only  small  amounts 
of  this  oxidizer  are  used  In  any  given  composition  due  to  the  constraint  H2/HCI  = 
97/3.  Best  (H,  + HCl)  purity  values  are  obtained  by  use  of  either  WO^  or  H2W0^ 
with  hydrogen-containing  fuels,  coolants  and  chlorides.  Both  of  these  oxidizers 
yield  the  metal  W(c)  (condensed  phase)  as  a combustion  product  which  can  be 
filtered  out  or  retained  In  the  motor,  to  prevent  contamination  of  the  working 
fluid.  Tungsten  Is  thermod3niamlcally  stable  at  the  temperatures  of  Interest 
(2800°K)  with  respect  to  H2 , HCl,  and  metal  oxides,  so  that  contamination  of  the 
working  fluid  with  tungsten  reaction  products  Is  not  a problem.  Other  possible 
oxidizers  Include  aluminum  nitrate,  A1(N02)2*9H20,  zirconium  nitrate,  Zr(N02)2‘ 
5H2O  and  zlrconyl  nitrate,  Zr0(N02)2 *2N20.  The  best  of  these  Is  Al-nltrate, 
yielding  2.4  mols  H2/IOO  gm,  which  Is  seen  to  be  Inferior  to  both  N2H^N02  and 
NH^NO^.  The  use  of  hydrates  as  propellant  Ingredients  Is  not  good  practice. 

The  organic  oxidizers,  NC  and  TMETN,  received  only  minor  attention  here  because 
contamination  of  the  H2  yield  with  CO  and  N2  Is  significant.  Calcium  and  mag- 
nesium nitrates  were  not  considered  because  of  the  Instability  at  2800*’K  of 
CaO(c)  and  MgO(d),  which  are  major  combustion  products  of  these  oxidizers. 

Both  oxides  react  significantly  with  H2  and  HCl  at  high  temperatures  to  yield 
H2O  and  other  contaminants  (Ca,  Mg,  CaOH,  MgOH,  CaCl,  CaCl2»  MgCl2).  Similarly, 
the  other  alkaline  earth  nitrates,  Sr(N02)2  and  Ba(lK)2)2»  which  are  well-known 
oxidizers  used  In  colored  flare  signals,  are  not  considered  usable  here.  Alkali 
metal  nitrates  and  perchlorates  (e.g.,  NaNO^,  CaClO^,  KNO^,  etc.)  yield  the 
chlorides  quantitatively  In  combustion  products.  These  species,  although  gaseous 
at  stagnation  or  chamber  conditions,  would  condense  In  the  cavity  and  for  this 
reason  these  oxidizers  are  not  considered  as  propellant  Ingredient  candidates. 
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3.2.2 


Fuels 


As  mentioned  previously,  hydrocarbons  are  not  desirable  as  fuels, 
because  combustion  to  CO  is  not  energetic  enough  to  yield  a flame  temperature 
of  2800®K,  and  combustion  to  CO2  yields  H2O  due  to  the  water-gas  reaction. 

For  the  attainment  of  high  temperature  without  H2®  contamination,  metallic 
fuels  must  be  employed. 

In  terms  of  H2  yield,  AlH^  > ZrH2  > A1  > Zr;  in  terms  of  H2  yield 
purity,  AIH^  > ZrH2  > A1  = 2r,  and  in  terms  of  heating  value,  A1  > Zr  > ZIH^  > 
ZrH2.  Since  AlH^  has  stability  and  sensitivity  problems,  it  is  not  a 
viable  candidate.  Consequently,  the  value  of  the  fuels  in  propellant 
compositions  la  considered  primarily  as  the  source  of  energy  rather  than 
of  H2«  Either  ZrH2  or  Zr  are  favored  in  preference  to  Al.  This  is  because 
Zr02,  the  major  Zr  combustion  product  at  the  OMOX  point,  is  thermodynamically 
more  stable  than  Al20^.  Even  though  the  use  of  Al  compares  favorably  with  Zr 
and  ZrH2  in  terms  of  H2  yield  purity,  A.120^  dissociates  somewhat  and  reacts 
with  H2O  as  a combustion  product.  In  Zr-contalnlng  propellants,  the  H2/H2O 
molar  ratio  is  'V250,  but  is  only  'V'80  in  the  combustion  products  of  propellants 
containing  Al.  The  additional  effect  of  this  phenomenon  on  the  combustion 
products  composition  at  throat  conditions  is  discussed  in  Section  3.3. 

For  reasons  stated  previously,  due  to  CaO  and  MgO  stability,  neither 
CaH2  nor  MgH2  were  considered  as  fuels. 

3.2.3  Coolants 

The  combustion  of  stoichiometric  mixtures  of  candidate  oxidizers 
and  fuels  yields  values  of  flame  temperature  greater  than  2800°K.  In  order  to 
reduce  the  flame  temperature  to  2800®K,  coolants  are  required  and  those  of 
choice  include  Zr(OH)^,  AlCOH)^,  (0112)2  HDB.  Stoichiometric  mixtures  of 
the  hydroxides  with  candidate  fuels  serve  as  coolants  to  obtain  a flame  tem- 
perature of  2800®K,  and  yield  H2  as  the  only  gaseous  combustion  product.  The 
as-received  form  of  Zr(OH)^  is  the  hydrate  Zr  (0H)^*nH_20  (n  = 4 to  5)  which 
must  be  dried  carefully  prior  to  use  (48  hours  at  145 ®F  in  a forced-draft 
oven  has  been  found  to  be  satisfactory).  Care  must  be  used  because  excessive 
drying  at  too  high  a temperature  drives  the  as-received  product  to  Zr02. 


- 47  - 


Stoichiometric  mixtures  of  (CH^)!!  and  candidate  oxldlzets  are  coolants 

If  the  oxidation  ratio  (0„)  of  the  mixture,  defined  as  the  molar  ratio  of 

K 

oxygen  to  carbon,  Is  unity.  At  values  of  0^^  <1,  soot  and  hydrocarbon  (HC) 
fragments  (CH,  CH2,  CH^,  CH^,  etc.)  appear  as  combustion  products,  while  at 
>1,  CO2  and  H2O  appear  as  combustion  products.  Note  that  while  the  H2 
yield  of  (CH2)  n is  second  only  to  HDB  In  terms  of  mols/100  gm,  the  molar 
purity  of  the  H2  yield,  assuming  C -»■  CO  is  only  50%,  (CH2  CO  + H2,  H2/ (H2  + c6) 
= 1/2).  Because  of  the  high  degree  of  purity  required, the  use  of  hydrocarbons  in 
candidate  propellants  was  given  only  minor  consideration.  The  purity  effect 
can  be  minimized  by  using  the  hydrocarbon  In  propellants  containing  no  oxygen 

■I 

for  burning  the  carbon.  However,  soot  and  HC  fragments  would  result.  It  Is 

possible  to  filter  out  the  soot,  depending  on  the  ability  of  the  composition  to 
form  a clinker  on  combustion  (If  the  composition  does  not  clink,  the  soot  will 
be  generated  as  angstrom-sized  particles  and  aggomerates  thereof,  which  will 
probably  be  very  difficult  to  remove).  Due  to  this  and  the  unknown  effect  that 
HC  fragments  would  have  on  gain  and  available  power,  hydrocarbons  were  ruled  out. 

Of  the  three  coolants,  only  HDB  Is  not  commercially  available  in 
large  quantities.  However,  synthesis  of  this  compound  is  straight  forward 
and  results  In  a product  of  good  yield  and  high  purity.  It  is  the  best  H2  gas 
generator  Ingredient  of  all  the  candidate  Ingredients.  It  is  capable  of  self- 
sustained  combustion,  with  an  adiabatic  flame  temperature  of  1877°K  (B2N2H2  -»■ 

5H2  + 2BN(c)).  The  compound  la  not  a good  candidate  for  systems  of  flame 
temperature  equal  to  3000®K  because  of  significant  dissociation  of  BN(c)  at 
this  temperature,  and  subsequent  recombination  in  the  throat,  resulting  In  a 
dirty  dust-laden  working  fluid.  However,  BN(c)  Is  essentially  imdissoclated 
at  2800*“K  and  the  H2  yield  Is  essentially  clean.  HDB  Is  a fairly  reactive 
solid,  and  has  exhibited  chemical  Incompatibility  with  many  propellant  Ingre- 
dients of  Interest.  Storage  and  use  la  limited  to  very  dry  facilities  (<40% 
relative  humidity  at  70®F) . Its  use  with  NH^CIO^  may  result  In  explosions  on 
Ignition,  though  this  seems  to  be  a phenomenon  dependent  on  HDB  purity  and  the 
absolute  and  relative  contents  of  each  component  In  the  propellant  composition. 

The  oxidizers  and  the  coolants  are  the  primary  H2  sources  In  the 
propellants.  The  requirement  (see  Section  3.1) that  (H2  + HCl)  yield  be  >5% 
by  weight  @ (H2/HCI)  » 97/3  (molar)  dictates  that  the  H2  yield  values  as 
tabulated  In  Table  3.1)  be  >1.59  mol/100  gm  of  propellant.  Disregarding  the 
organics  and'AlH^,  It  Is  seen  that  the  only  oxldlzer-fuel  combinations  that 
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qualify  are  N2H^NQj  + either ZrH2  or  Al.  The  only  coolant  that  qualifies  is  HDB. 
Of  the  two  sources,  HDB  is  obviously  superior,  not  only  in  terms  of  H2 
yield,  but  also  H2  gas  purity. 

3.2.4  Chlorides 

The  chlorides  are  Included  in  propellants  as  Cl-donors  to  yield  HCl 
in  the  combustion  products.  All  are  commercially  available  and  easily  used 
and  handled.  The  H-contalnlng  compounds  also  yield  N2  as  a combustion  product, 
but  this  is  not  too  Important  because  the  chloride  content  of  propellants  is 
not  high  due  to  the  requirement  for  a high  value  of  H2/HCI  (97/3).  For 
compositions  in  which  no  contaminants  are  permitted  in  the  combustion  products, 
ZrCl^  can  be  used  as  the  Cl-donor.  For  example,  the  equimolar  mixture  of 
ZrCl^  and  Zr(OH)^  yields  HCl  and  Zr(^(c)  as  the  only  combustion  products. 

In  terms  of  H2  yield,  N2H5CI  > NH^Cl  > (NH^CIO^  + ZrH2)  > (NH^CIO^  + Al) 
> (NH^CIO,^  + Zr).  In  terms  of  HCl  yield,  NH^Cl  > N2H^C1  > (NH^CIO^  + Al)  > 

(NH^  CIO^  + Zr)  = (NH^  CIO^  + ZrH2).  Values  of  (H2  + HCl)  yield  purity  range 
from  74%  for  N2H2CI  up  to  88%  for  NH^CIO^  + ZrH2.  The  NH^ClO^-fuel  combinations 
are  of  course  more  energetic  than  either  of  the  chlorides. 
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3. 3 Calculation  of  Theoretical  Yield 

3.3.1  Procedure 

The  selection  of  candidate  propellant  compositions  Is  based, primarily 
on  computed  values  of  theoretical  flame  temperature  and  desired  combustion 
product  composition.  Basically  the  propellant  can  be  constructed  from  any 
combination  of  oxidizer,  fuel,  coolant  and  chloride  from  the  binary  combina- 
tions In  Table  3.1  by  means  of  two  equations  for  a heat  and  molar  ratio 
balance;  these  equations  are: 

rCXAH^^.go^l  - 0 Heat  Balance  (D 

E(XHCl)j,  HCl  Molar  Ratio  Balance  (2) 

where  are  multipliers.  In  (1)  the  condition  Is  that  combustion  is  adia- 
batic, while  In  (2),  H2/HCI  - 97/3  generally.  In  all  cases,  there  are  three 
multipliers  (e.g.,  one  for  oxldlzer-fuel,  one  for  fuel-coolant,  and  one  for 
chloride)  any  one  of  which  Is  arbitrary  and  can  be  set  to  unity.  Simultaneous 
solution  of  (1)  and  (2)  for  the  other  two  multipliers  and  normalization  to 
100  result  In  propellant  compositions  that  yield  the  flame  temperature  and  H2 
and  HCl  very  close  to  the  desired  values.  Equilibrium  calculations  of  resulting 
compositions  yield  slightly  different  values  than  predicted  by  using  (1)  and  (2) 
because  of  mixing  effects  on  entropy  and  the  resultant  equilibria.  If  required, 
compositions  are  then  sometimes  varied  slightly  to  attain  the  desired  results. 
For  example,  consider  the  system  N2H^N0g  + Zr  H2  + Zr(OH)^  + NH^CIO^.  What 
is  the  OMOX  composition  with  a flame  temperature  of  2800®K  and  H2/HCI  = 97/3? 
Consider  the  three  binaries: 

A.  40.467%  N2H^N03  + 59.533%  Zr  H2,  4H“ggQ  = -40.91, 

H2  " 1.685,  HCl  = 0 

B.  63.078%  Zr(OH)^  + 36.922%  Zr  H2,  AHjggg  = 30.48, 

H2  - 1.181,  HCl  - 0 

C.  38.654%  NH^CIO^  + 61.346%  ZrH2,  AHjggg  - -67.44, 

H2  » 1.144,  HCl  - .3048 
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and  set  x,  = 1 
A 


Let  X^,  Xg  and  X^  be  the  multipliers  for  compositions  A,  B and  C,  respectively, 
A=  1- 
From  (1) 

- 40.91  + 30.48X_  - 67.44X^  »=  0 

o L 


From  (2) 


1.685  + l.lSlXg  + 1.144Xg 

0.3048X^  “ “3 


From  vdilch  we  determine 


Xg  = 2.5287 
X^  = 0.53625 


The  resultant  composition  Is  comprised  of  lOOX^  ■ 100  gm  of  Composition  A plus 

lOOXg  = 252.87  gm  of  Composition  B and  lOOX^  >=  53.625  gm  of  Composition  C,  and  the 

weight  percentages  of  each  Ingredient  are  obtained  by  normalizing  to  the  sum 

X.  + X„  + X„  = 4.0650: 

ABC 

N2H^N03  - 40.467X^/(X^  ^C^ 

- 40.467/4.0650 


“ 9.955 

ZrH2  - (59.533X^  + 36.922Xg  + 61.346X(,)/ (X^  + Xg  + X^^) 
= 185.79  / 4.065 
= 45.706 

Zr(OH)^  - 63.078Xg/(X^  + Xg  + X^,) 

“ 159.50/4.065 
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39.239 


NH^CIO^  - 38.654Xj,/(X^  + Xg  + X^) 
- 20.728/4.065 


» 5.100 


Slight  variations  In  the  composition  to  compensate  for  mixing  effects,  as 
above,  are  generally  derived  by  solving  for  NH^CIO^  content  ([NH^CIO^]) 
around  the  above  value.  Thus,  replace  (2)  with 


38.654X^ 

x^TxT+x^ 


[NH^CIO^] 


and  solve  for  selected  values  of  [NH^CIO^],  say  [NH^CIO^]  = 3,  4,  5,  6 and 
7%. 

Calculations  of  theoretical  flame  temperature  and  combustion  product 
composition  were  carried  out  for  adiabatic  combustion  at  750  psl  followed  by 
Isentroplc  expansion  of  the  gaseous  combustion  products  to  throat  conditions. 

The  condensed  species  (^20^,  Zr02,  W,  BN,  etc.)  In  these  calculations  were 
assumed  to  be  retained  In  the  combustion  chamber  and  not  allowed  to  undergo 
the  expansion  process  with  the  gases.  Tabulated  values  for  H2  yield,  HCl 
yield,  and  other  gaseous  species  are  throat  values. 

The  results  of  these  calculations  and  the  finalized  compositions, 
for  H2/HCI  ^ 97/3  are  summarized  In  Tables  3.3  to  3.10.  Results  are  discussed 
below. 

3.3.2  Discussion  of  Results 

The  results  of  the  thermochemistry  calculations  are  summarized  by 
oxidizer  for  compositions  containing  NH^NO^  (Tables  3.3  and  3.4),  N2H2NO2 
(Tables  3.5  and  3.6),  NH^CIO^  (Table  3.7),  H2W0^  (Table  3.8),  WO^  (Table  3.9) 
and  Zr0(N0g)2 • 2H2O  (Table  3.10).  Some  calculations  were  carried  out  for  temper- 
atures other  than  2800‘’K  and  these  also  are  presented  In  the  Tables,  but  without 
discussion. 
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3 . 3 . 2 . 1 Compositions  Containing  NH^NO 


The  results  for  compositions  containing  either  Zr  or  ZrH2  are  shown 
In  Table  3.3  whereas  results  for  either  A1  or  AlH^  are  shown  In  Table  3.4. 

All  of  the  compositions  containing  organic  coolants  (CH2  or  HTPB)  easily  meet 
the  (H2  + HCl)  weight  yield  requirement  and  are  easily  reduced  to  practice, 
but  are  very  poor  from  the  standpoint  of  (H2  + HCl)  gas  purity.  The  composi- 
tion that  contains  off-the-shelf  Ingredients  and  that  comes  closest  to  meeting 
the  yield  requirement  (>90%)  Is  3.3.4,  which  Is  comprised  of  ZrH2, 

Zr(OH)^,  and  NH^CIO^.  This  composition  was  selected  as  one  of  the  candidates 
and  is  considered  to  be  representative  of  the  best  that  can  be  done.  In  terms 
of  yield  and  purity,  with  off-the-shelf  ingredients.  Notice  though  that  the 
use  of  HDB  as  a coolant  (Compositions  3.3.21  arid  3,3.22)  greatly  Increases 
both  yield  and  gas  purity. 

One  of  the  most  serious  drawbacks  to  the  use  of  A1  or  AlH^  as  a fuel 
is  the  recombination  of  reaction  products  in  the  throat.  As  noted  previously, 
AI2O2  Is  less  stable  than  Zr02  at  the  desired  stagnation  conditions  which  re- 
sult in  the  reactions 

^1303(0)  + 2H2  + 2HC1  2A1C1  + 3H2O  (1) 

^1203(0)  + H2  + 4HC1  -»■  2AICI2  ^”2° 

In  the  combustion  chamber.  The  reverse  of  these  reactions  occurs  at  the  throat. 
Although  only  gaseous  species  are  expanded  to  throat  conditions,  recombination 
of  AlCl  or  AICI2  with  H2O  results  in  small  amounts  of  Al202(c)  formation  at 
throat  conditions.  The  yield  of  AI2O2  through  recombination  Is  1-2%  by 
weight  of  the  working  fluid.  Since  the  nozzle  does  not  allow  time  for  equili- 
brium to  be  obtained,  these  are  strictly  thermodynamic  results  and  the  amount 
of  AI2O2  that  actually  forms  and  the  extent  of  interference  this  phenomenon 
may  have  on  nozzle  and/or  laser  performance  is  not  known.  A similar  phenomenon 
Is  found  for  HDB-contalnlng  propellants  due  to  recombination  of  the  decomposi- 
tion and  reaction  products  of  BN(c).  However,  this  compound  can  be  used  In  a 
two-compartment  motor  system,  and  the  BN  can  be  filtered  prior  to  admixture  of 
the  H2  with  the  other  combustion  products.  In  this  way,  the  recombination 
phenomenon  Is  obviated. 
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Table  3.3 

Theoretical  Values  of  Flame  Temperature  and  Yield  for  Selected  OMOX  Compositions: 
Compositions  Containing  NH^NO^  and  Either  Zr  or  ZrH^ 


Yield 


Comp 

Propellant  Composition,  Wt 

% 

T “K  * 

f’ 

H2/HCI* 

(H^+HCl) * 

ID 

Oxidizer 

Fuel 

Coolant 

Chloride 

0750  psia 

(molar) 

Wt  % 

Vol  % 

3.3.1 

29.144 

NH.NO- 

57.325 

ZrH- 

10.948 

Zr(OH), 

2.583 

NH.Cl 

3000 

32.1 

4.85 

79.4 

3.3.2 

17.950 

4 J 

49.203 

z 

30.482 

h 

2.365 

4 

2799 

32.2 

4.46 

85.1 

3.3.3 

20.241 

NH.NO- 

54.914 

ZrH_ 

19.845 

Zr(OH), 

5 

NH  CIO. 

3000 

33.1 

4.35 

83.0 

3.3.4 

10.482 

4 j 

47.132 

L 

38.386 

4 

2801 

38.2 

3.84 

89.5 

3.3.5 

4.008 

42.838 

48.440 

4.714 

2400 

30.7 

3.93 

94.6 

3.3.6 

36.955 

NH.NO- 

58.194 

ZrH- 

1.917 

(CH  ) 

2.934 

NH.Cl 

3000 

32.2 

5.51 

73.5 

3.3.7 

39.538 

4 J 

52.007 

L 

5.137 

z n 

3.318 

2801 

32.3 

6.27 

69.4 

3.3.8 

40.574 

49.532 

6.426 

3.468 

2602 

32.3 

6.57 

68.2 

3.3.9 

34.015. 

NH,N0_ 

56.126 

ZrH- 

3.859 

(CH  ) 

6 

NH.CIO 

3000 

36.3 

6.45 

70.9 

3.3.10 

35.990 

49.702 

L 

7.308 

Z u 

7 

2802 

34.2 

6.26 

67.3 

3.3.11 

37.162 

47.145 

8.693 

7 

2602 

35.9 

6.47 

66.2 

3.3.12 

33.918 

NH.NO- 

54.736 

ZrH- 

4.655 

HTPB 

6.691 

NH  CIO 

3000 

31.0 

5.58 

69.0 

3,3.13 

36,803 

H J 

47,670 

z 

8.527 

7 

H ■ H 

2802 

33.1 

6.11 

64.4 

3,3,14 

38.122 

44,740 

10.137 

7 

2602 

34.5 

6.30 

63.0 

3.3.15 

39,461 

41,768 

11.771 

7 

2402 

35.9 

6.47 

. 61.5 

3.3.16 

39,628 

NH.NO- 

43.982 

Zr 

10.394 

HTPB 

6 

NH  CIO 

3000 

31.7 

5.07 

56.1 

3.3.17 

41.904 

4 j 

38,876 

13.220 

6 

2801 

35.6 

5.51 

55.2 

3.3.18 

42.118 

36,226 

14,656 

7 

2601 

32.1 

6.02 

54.9 

3.3.19 

43.208 

33.784 

16,008 

7 

2401 

33.7 

6.21 

54.5 

3.3.20 

23.507 

NH,N0_ 
4 3 

40,182 

Zr 

30,311 

HDB 

6 

NH.Cl 

4 

2825 

34.0 

10.1 

89.3 

3.3.21 

2.482 

NH4NO3 

33.744 

Zr 

44.774 

HDB 

19 

NH^CIO^ 

2812 

32.0 

12.7 

95.3 

3.3.22 

7.782 

NH.NO. 
4 3 

34.896 

Zr 

+ 

40,322 

1 

HDB 

16 

NH.CIO, 
4 4 

2818 

34.1 

11.8 

92.6 

* Notes; 

Tj  is 

adiabatic 

flame  temperature  at  a combustion 

pressure 

of  750  psia. 

H2/HCI  and  yield  (H^+HCl)  refer  to  combustion  product  values  at  throat  conditions. 
Yield  (H2+HCI) , vol.  % is  the  yield  purity  value. 


Table  3.4 

Theoretical  Values  of  Flame  Temperature  and  Yield  for  Selected  OMCK  Compositions: 
Compositions  Containing  NH,N0_  and  Either  Al  or  A1H„. 


Comp 

Propellant  Composition,  Wt 

% 

ID 

Oxidizer 

Fuel 

Coolant 

3.4.1 

23.140 

NH,N0_ 

30.571 

Al 

43.289 

Al(OH), 

3.4.2 

18.016 

29.462 

50.069 

3.4.3 

13.647 

NH  NO 

30.219 

Al 

50.119 

Al(OH) , 

3.4.4 

10.212 

H J 

29.179 

55.609 

3.4.5 

3.4.6 

59.033 

58.973 

NH,N0„ 
4 3 

26.953 

25.097 

Al 

10.014 

11.430 

(CH2>^ 

3.4.7 

59.442 

23.294 

13.082 

3.4.8 

3.4.9 

52.190 

52.463 

NH^N03 

25.431 

23.541 

Al 

12.379 

13.996 

3.4il0 

53.005 

21.791 

15.502 

3.4.11 

52.234 

NH,N0- 

23.772 

Al 

13.994 

HTPB 

3.4.12 

53.454 

21.870 

15.676 

3.4.13 

53.707 

19.947 

17.346 

3.4.14 

53.952 

18.100 

18.948 

3.4.15 

36.026 

NH,N0_ 
4 3 

35.610 

AIH3 

22.364 

Al(OH), 

3.4.16 

39 . 289 

NH  NO 

35.162 

AlH- 

18.549 

Al(OH). 

+ 1 

CH2 

Yield 


Chloride 

Tf.  °K, 
@750  psia 

H^/HCl 

(molar) 

(H2+HCI) 

Wt  % Vol  % 

3 

NH.Cl 

2989 

33.8 

4.56 

80.6 

2.453 

4 

2791 

37.8 

4.39 

84.8 

6.015 

NH.CIO, 

2995 

34.4 

4.19 

85.3 

5 

4 4 

2794 

38.2 

4.10 

89.2 

4 

NH.Cl 

2971 

39.4 

6.67 

59.9 

4.5 

4 

2776 

33.0 

7.48 

60.2 

4.182 

2584 

34.6 

7.73 

60.0 

10 

NH  CIO 

2970 

34.4 

7.02 

58.6 

10 

4 4 

2777 

32.9 

7.54 

59.0 

9.702 

2585 

33.0 

7.92 

59.0 

10 

NH.CIO, 

2970 

32.2 

6.76 

55.5 

9 

4 4 

2778 

34.3 

6.97 

55.5 

9 

2585 

33.2 

7.36 

55.2 

9 

2393 

32.6 

7.68 

54.8 

6 

NH.Cl 

4 

2784 

33.9 

10.1 

85.7 

6 

NH.Cl 

2783 

34.5 

10.2 

83.6 

3 . 3 . 2 . 2 Compositions  Containing  N^H^NOj 

These  results  are  shown  In  Table  3.5  for  compositions  containing 
either  Zr  or  ZrH2  as  a fuel  and  In  Table  3.6  for  A1  and  AlH^  as  fuels.  Com- 
pared to  NH^NO^,  the  use  of  this  oxidizer  In  comparable  formulations  results 
In  slight  Increases  In  yield  and  decreases  In  gas  purity.  For  example,  compare 

3.3.2  and  3.5.2,  at  2800 ®K  the  composition  containing  NH^  has  a weight  yield 
of  4.46%  and  a purity  of  85.1%,  idiereas  the  composition  containing  N^H^NO^ 
has  a weight  yield  of  4.64%  and  a purity  of  81.4%. 


3. 3. 2. 3 Compositions  Containing  Either  H^WO^^  or  WO^ 

These  oxidizers  are  Ideal  for  use  In  thermite  systems,  resulting  In 
excellent  clinkers.  Both  result  In  high  values  of  (H2  + HCl)  gas  purity  since 
no  N2  Is  Introduced  In  the  combustion  products.  When  used  In  propellants 
made  with  off-the-shelf  Ingredients,  H2W0^  la  preferred.  The  composition 
3.8.1,  comprised  of  H2W0^,ZrH2,  Zr,  and  WH^CIO^  Is  representative  of  the  best 
that  can  be  attained  from  thermite  systems  made  with  off-the-shelf  materials 
and  was  selected  as  one  of  the  propellant  candidates.  Note  that  the  yield  of 
this  composition,  2.59%,  Is  very  low;  however.  It  Is  also  very  pure.  When 
used  with  HDB,  there  Is  not  much  advantage  between  the  two  oxidizers  and.  In 
fact,  the  performance  Is  comparable  to  the  use  of  NH^NO^  or  N2H^N02  (>10% 
weight  yield).  The  composition  selected  as  one  of  the  candidates  was  composi- 
tion 3.8.8. 
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Table  3.5 

Theoretical  Values  of  Flame  Temperature  and  Yield  for  Selected  OMOX  Compositions: 
Compositions  Containing  N2H^N02  and  Either. Zr  or  ZrH2 


Propellant  Composition,  Wt  % 


T^,  “K,  H2/HCI 


Yield 

(H2+HCI) 


ID 

Oxidizer 

Fuel 

Coolant 

Chloride 

@750  psia 

(molar) 

Wt  % 

Vol 

3.5.1 

31.997 

N2H3NO3 

53.994 

ZrH2 

11.484 

ZrCOH), 

4 

2.725 

NH.Cl 

4 

3000 

32.1 

5.11 

75*1 

3.5.2 

19.832 

47.097 

30.615 

2.456 

2800 

32.2 

4.64 

81.4 

3.5.3 

31.621 

N2H3NO3 

53.292 

ZrH2 

11.571 

Zr(OH)^ 

3.516 

N2H3CI 

3000 

32.1 

5.15 

74.6 

3.5.4 

22.126 

N2H5NO3 

52.444 

ZrH2 

20.430 

Zr(OH)^ 

5 

NH^CIO^ 

3000 

34.4 

4.46 

79.2 

3.5.5 

41.560 

N2H3NO3 

52.696 

ZrH2 

2.541 

^^2>n 

3.203 

NH^Cl 

2950 

32.2 

6.02 

68.9 

3.5.6 

43;917 

N2H3N03 

47.368 

5.188 

3.527 

2801 

32.3 

6.66 

66.3 

3.5.7 

40.759 

N2H5NO3 

53.131 

ZrH2 

2.056 

4.054 

N2H3CI 

3000 

32.2 

5.94 

68.9 

3.5.8 

37.070 

N2H5NO3 

51.781 

ZrH2 

4.164 

6.985 

NH^CIO^ 

3000 

32.2 

5.97 

67.5 

3.5.9 

35.114 

N2H3NO3 

51.658 

ZrH2 

9.290 

HDB 

3.938 

NH^Cl 

2999 

32.1 

7.38 

79.4 

3.5.10 

25.758 

N2H3NO3 

37.074 

Zr 

29.168 

HDB 

8 

N2H3CI 

2826 

32.7 

10.4 

85.7 

3.5.11 

3.969 

N2H3NO3 

33.661 

Zr 

44.370 

HDB 

18 

NH^CIO^ 

2817 

33.2 

12.6 

94.8 

3.5.12 

23.808 

N„H-N0- 

34.267 

Zr 

31.825 

HDB 

6 

ZrCl 

2831 

36.8 

9.84 

88.7 

+4.100 

Zr(OH) 

4 

Theore 


Comp  Propellant  Composition,  Wt  % 

ID  Oxidizer  Fuel Cool. 


3.6 

1 

25.786 

N2H5N03 

29.177 

A1 

42.037 

3.6 

2 

20.438 

28.294 

48.268 

3.6 

3 

25.523 

N2H5NO3 

28.876 

A1 

41.601 

3.6 

4 

14.546 

N2H5NO3 

29.383 

A1 

49.571 

3.6 

5 

62.206 

N2H5NO3 

23.884 

A1 

8.910 

3.6 

6 

63.059 

21.867 

li.861 

3.6 

7 

61.228 

N2H5NO2 

23.477 

A1 

8.796 

3.6 

8 

54.611 

N2H5NO3 

22.537 

A1 

11.852 

3.6 

9 

36.180 

N2H3NO3 

20.536 

A1 

34.784 

3.6 

10 

47.622 

N2H3NO3 

34.316 

AIH2 

11.062 

3.6 

11 

37.348 

N2H3NO3 

32.764 

AIH3 

23.888 

3.6 

12 

24.668 

N2H5NO3 

30.736 

AIH3 

39.588 

3.6 

13 

57.759 

N2H3NO3 

32.379 

Alfl3 

2.862 

3.6 

14 

59.143 

N2H5NO3 

29.044 

AIH3 

5.813 

3.6 

15 

18.264 

N2H3NO3 

33.600 

AIH3 

36.136 

3.6 

16 

6.248 

N2H3NO3 

31.246 

AIH3 

51.506 

3.6 

17 

32.337 

N2H3NO3 

28.895 

AIH3 

33.768 

3.6 

18 

12.474 

N2H3NO3 

27.300 

AIH3 

50.226 

JL 


Tf,  “K, 

H^/HCl 

(H2+HCI) 

mt 

Chloride 

@750  psia 

(molar) 

Wt  % 

Vol  % 

AKOH)^ 

3 

NH.Cl 

4 

2989 

35.8 

4.70 

76.6 

3 

2789^ 

32.0 

4.86 

80. 8 

AKOH)^ 

4 

N,H,.C1 

2988 

34.6 

4.81 

75.8 

AKOH)^ 

6.5 

NH.CIO, 
4 4 

2995 

32.5 

4.39 

82.6 

5 

KH.Cl 

4 

2969 

32.4 

7.47 

58.4 

4.213 

2777 

37.4 

7.57 

58.4 

6.5 

N2H^C1 

2968 

32.1 

7.54 

■58.3 

11 

NH4CIO4 

2967 

32.5 

7.49 

57.3 

HDB 

8.5 

NH^Cl 

2950 

32.0 

12.8 

84.5 

AKOH)^ 

7 

NH.Cl 

4 

2982 

32.6 

10.44 

78.7 

AKOH)^ 

6 

NH^Cl 

2746 

33.1 

10.1 

82.8 

AKOH)^ 

5.108 

NH,  Cl 
4 

2391 

34.0 

9.52 

87.7 

AKOH)^ 

7 

NH^Cl 

2978 

34.2 

10.8 

73.5 

6 

NH^Cl 

2780 

37.2 

10.7 

71.5 

akoh)^ 

12 

NH.CIO, 
4 4 

2794 

34.3 

9.20 

88.8 

AKOH)^ 

11 

NH,C10, 
4 4 

2395 

32.2 

8.99 

95.0 

Zr(OH), 

4 

5 

NH^Cl 

2784 

35.6 

8.63 

82.9 

Zr(OH), 

4 

10 

NH,C10, 
4 4 

2794 

33.5 

7.56 

90.0 

[ 


! 


Table  3.7 

Theoretical  Values  of  Flame  Temperature  and  Yield  for  Selected  OMOX  Compositions: 


Miscellaneous  Compositions  Containing  NH^CIO^ 

Comp  Propellant  Composition,  Wt  % 

ID  Oxidizer  Fuel  Coolant  Chloride 


3.7.1 

42.794 

NH.CIO, 
4 4 

26.206 

A1 

3.7.2 

16.725 

NH.CIO, 
4 4 

9.717 

A1 

3.7.3 

13.317 

NH,C10, 
4 4 

20.677 

Zr 

3.7.4 

13.586 

NH,C10, 
4 4 

19.921 

Zr 

3.7.5 

3.095 

NH,C10, 
4 4 

37.850 

Zr 

3.7.6 

3.000 

NH,C10, 
4 4 

35 . 728 
+2.596 

Zr 

ZrH. 

3.7.7 

16.714 

NH, CIO, 
4 4 

26.526 

ZrH, 

31 

HDB 

70.161 

HDB 

3.397 

ZrCl 

66.006 

HDB 

63.493 

HDB  . 

3 

ZrCl 

0.205 

(CH2)^ 

58.676 

Zr(OH)^ 

56.760 

HDB 

Yield 


Tf,  “K, 

H2/HCI 

(H  +HC1) 

0750  psia 

(molar) 

Wt  Z 

Vol  : 

2798 

42.4 

9.28 

76.5 

2347 

34.4 

18.6 

98.2 

2394 

50.6 

15.6 

98.8 

2393 

32.9 

17.0 

98.7 

2800 

30.2 

2.55 

96.2 

2800 

31.6 

2.5i 

98.0 

2384 

37.1 

15.7 

98.4 

I 

Ln 

VO 


i 


i 


I 


r 


f 


o\ 

o 


Table  3.8 

Theoretical  Values  of  Flame  Temperature  and  Yield  for  Selected  OMOX  Compositions: 
Compositions  Containing  H2W0^  and  Either  Zr  or  ZrH2 


Propellant  Composition,  Wt  % 


‘f’ 


’K. 


H2/HCI 


Yield 

(H2+HCI) 


ID 

Oxidizer 

Fuel 

Coolant 

Chloride 

@750  psia 

(molar) 

Wt  % 

Vol 

3.8.1 

40.879 

HzWO^ 

42.968 

ZrH2 

13.153 

Zr(OH)^ 

3 

NH4CIO4 

2799 

32.6 

2.59 

98.0 

3.8.2 

20.641 

H2WO4 

41.453 

ZrH2 

34.051 

Zr(OH)^ 

3.855 

NH^CIO^ 

2400 

30.8 

3.21 

98.3 

3.8.3 

54.657 

H2WO4 

42.018 

ZrH2 

0.825 

(CHp^ 

2.5 

NH,  CIO . 
4 4 

2800 

35.9 

2.30 

91.4 

3.8.4 

57.018 

H2WO4 

41.649 

ZrH2 

0.127 

<“2>n 

1.206 

ZrCl, 

4 

2798 

33.5 

2.09 

98.2 

3.8.5 

55.282 

H2W0^ 

41.856 

ZrH2 

0.615 

0.465 

+1.782 

ZrCl, 

NH4CIO4 

2799 

31.9 

2.31 

93.1 

3.8.6 

32.348 

H2WO4 

30.354 

4^.827 

ZrH 

Zr 

11.363 

+9.292 

HDB 

Zr(OH)^ 

6.816 

NH4CIO4 

2780 

34.6 

4.92 

96.4 

3.8.7 

47.996 

H2WO4 

35.044 

Zr 

14.533 

HDB 

2.427 

NH^Cl 

2802 

41.5 

4.25 

96.5 

3.8.8 

11.821 

33.475 

Zr 

38.704 

HDB 

16 

NH.CIO. 
4 4 

2808 

33.3 

10.8 

95.9 

3.8.9 

21.952 

H2WO4 

33.623 

Zr 

30 

+1 

HDB 

13.425 

NH.CIO, 
4 4 

2802 

32.2 

8.91 

93.7 

3.8.10 

48.988 

H2WO4 

7.096 

+33.816 

Zr 

Zrfi2 

6.600 

+1 

Zr(OH), 

(®2>n 

2.5 

NH,  CIO, 
4 4 

2800 

35.2 

2.26 

89.9 

I 


Table  3.9 

Theoretical  Values  of  Flame  Temperature  and  Yield  for  Selected  OMOX  Compositions: 

Compositions  Containing  WO-  and  Either  Zr  or  ZrH- 

^ Yield 


Comp 

ID 

Oxidizer 

Propellant  Composition,  Wt  % 

Fuel  Coolant 

Chloride 

Tf,  °K, 
@750  psia 

H2/HCI 

(molar) 

(H2+HCI) 

Wt  % Vol  % 

3.9.1 

45.181 

WO3 

38.685 

ZrH^ 

14.124 

Zr(OH)^ 

2 

NH.CIO, 
4 4 

2789 

36.4 

1.86 

98.2 

3.9.2 

60.983 

WO3 

36.097 

ZrH^ 

.1.084 

(CH2)^ 

1.836 

NH.CIO, 
4 4 

2800 

31.4 

1.55 

85.1 

3.9.3 

62.957 

WO3 

35.697 

ZrH2 

0.596 

<®2>n 

0.750 

ZrCl, 

4 

2799 

33.0 

1.30 

90.6 

3.9.4 

62.143 

WO3 

11.593 

Zr 

25.637 

ZrH2 

0.627 

N2H3CI 

3000 

32.2 

0.916 

95.9 

3.9.5 

60.239 

WO3 

7.178 

Zr 

31.199 

ZrH2 

1.384 

NH^CIO^ 

3000 

30.2 

1.13 

97.2 

3.9.6 

9.960 

WO3 

28.959 

Zr 

49.327 

49.254 

Zr(OH) . 
ZrH2 

2.5 

NH^eiO^ 

2800 

35.4 

2.28 

98.1 

3.9.7 

53.819 

WO3 

3-1.763 

Zr 

12-845 

HDB 

1.846 

NH^Cl 

2998 

32.2 

3.45 

97.4 

3.9.8 

20.113 

WO3 

32.056 

Zr 

34.831 

HDB 

13 

NH4CIO4 

2811 

36.2 

9.27 

96.0 

3.9.9 

17.840 

WO3 

10.529 

+20.032 

Zr 

ZrH2 

37.889 

HDB 

13 

NH,C10, 
4 4 

2600 

34.2 

11.0 

98.1 

!-■ 

I 


Table  3.10 

V, 

;?  Theoretical  Values  of  Flame  Temperature  and  Yield  for  Selected  OMOX  Compositions: 

Compositions  Containing  ZrO(NO^)  •2H„0  and  Either  Zr  or  ZrH- 

J z z Z Yield 


Comp  Propellant  Composition,  Wt  % ^f’  H^/HCl-  (H2+HCI) 


ID 

Oxidizer 

Fuel 

Coolant 

Chloride 

@750  psia 

(molar) 

Wt  % 

Vol  % 

3.10.1 

34.593 

ZrO(N02)2*2H2O 

50.168 

ZrH2 

13.549 

Zr(OH), 

4 

1.690 

NH.Cl 

4 

3000 

32.1 

3.17 

86.7 

3.10.2 

21.348 

Zr0(N02)2*2H20 

44.812 

ZrH2 

32.024 

Zr(OH)  , 
4 

1.816 

NH.Cl 

4 

2799 

32.2 

3.43 

91.5 

3.10.3 

26.311 

Zr0(ND3)2'.2H20 

49.998 

Zr«2 

19.691 

Zr  (OH)^ 

4 

NH.CIO, 
4 4 

3000 

30.3 

3.29 

89.1 

3.10.4 

12.840 

Zr0(N02)2*2H20 

44.598 

ZrH2 

38.562 

Zr(OH)^ 

4 

NH.CIO, 
4 4 

2800 

32.6 

3.46 

94.0 

3.10.5 

46.846 

ZrO(N02)2*2H2O 

48.784 

ZrH2 

2.532 

1.838 

NH.Cl 

4 

3000 

32.2 

3.45 

74.4 

3.10.6 

50.069 

ZrO(N0^)2*2H2O 

42.045 

ZrH2 

5.743 

2.143 

NH.Cl 

4 

2802 

32.3 

4.05 

68.0 

3.10.7 

44.307 

Zr0(N02)2*2H20 

47.924 

ZrH2 

3.769 

4 

NH.CIO, 
4 4 

3000 

34.2 

3.56 

71.8 

3.10.8 

46.725 

ZrO(N03)2*2H20 

41.084 

ZrH2 

7.191 

5 

NH,C10, 
4 4 

2802 

32.1 

4.28 

66.1 

3.4 


Properties  of  Candidate  Propellants 


3.4.1  General 

The  three  candidate  compositions  chosen  for  evaluation  in  motor 
firings  are  Compositions  3.3.4,  3.8.1  and  3.8.8.  The  first  two  selections  are 
about  the  beat  that  can  be  obtained  with  off-the-shelf  Ingredients  with  regard 
to  yield  and  purity,  while  the  last  is  a composition  containing  HDB.  As  indi- 
cated in  the  tables  of  the  previous  section,  there  are  other  propellant  selec- 
tions besides  3.8.8  that  have  a high  yield;  however,  only  one  high-yield  high- 
risk  system  was  selected.  This  selection  was  based  on  the  system  that  appeared 
to  have  the  highest  probability  of  success.  The  Composition  3.3.4  is  comprised 
of  NH^NO^,  ZrH2,  Zr(OH)^,  and  NH^CIO^;  Composition  3.8.1  is  comprised  of  H2W0^, 
ZrH2,Zr(0H)^,  and  NH^ClO^,  and  Composition  3.8.8  is  comprised  of  H2W0^,  Zrj 
HDB  and  NH^CIO^.  The  safety  and  ballistic  properties  of  each  composition  were 
measured,  and  are  given  in  Table  3.11.  Detailed  theoretical  results  for  flame 
and  throat  temperatures  for  the  combustion  product  composition  at  a chamber 
pressure  of  750  psla  are  given  for  each  of  the  compositions  in  Table  3.12. 
Burning  rate  curves  of  rate  vs  pressure  are  shown  for  compositions  3.3.4  and 

3.8.1  in  Figures  3.2  and  3.3,  respectively. 

All  of  the  compositions  are  prepared  by  first  blending  the  ingredients 
and  then  pressing  them  to  the  desired  size  and  shape  at  30  KSI.  For  each 
composition,  two  values  of  the  nozzle  discharge  coefficient  are  given,  one  for 
the  gas-only  flow  condition  (which  is  the  design  condition)  and  the  other  for 
the  flow  of  both  gases  and  solids  (given  in  parenthesis) . Burning  rates  are 
for  specimens  ~0.5"  dia  x ~l"long  pressed  at  30  KSI  in  steel  tubes.  Burning 
rates  derived  from  subscale  motor  firings  (see  Section  3.5)  are  also  shown. 
Results  are  discussed  below. 

3.4.2  Properties  of  Composition  3.3.4 

This  composition  is  the  best  (in  terms  of  gas  purity)  of  the  nitrate- 
containing  systems  comprised  of  off-the-shelf  ingredients.  The  composition 
presses  up  fairly  well  at  30  KSI,  ignites  easily  and  bums  uniformly  in  the 
pressure  range  400-1500  psia.  Combustion  at  750  psia  yields  combustion  products 
comprised  of  8.04%  by  weight  gases  and  91.96%  solids.  At  this  pressure,  and 
at  higher  values,  the  solids  form  a good  clinker.  At  low  pressure  (<100  psia) 
the  condensed  combustion  products  do  not  adhere  and  are  entrained  by  the  gases. 
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Table  3.11 

PROPERTIES  OF  CANDIDATE  PROPELLANTS 


Property 

Comp  3.3.4 

Comp  3.8.1: 

Comp  3.8.8 

A. 

Composition,  Wt% 

• 

Oxidizer 

10.482  NH^NO^ 

40.879  H2W0^ 

11.821  H2W0^ 

Fuel 

47.132  ZrH2 

42.968  ZrH2 

33.475  Zr 

Coolant 

38.186  Zr(OH)^ 

13.153  Zr(Ott)^ 

38.704  HDB 

Chloride 

4.000  NH^CIO^ 

3.000  NH^CIO^ 

16.000  NH.CIO, 
4 4 

B. 

Theoreticals 

Flame  Tempera- 
ture, °K 

2799 

2799 

2808 

(H^/Mci)  molar 
Yield  (H2+HCI) 

38.7 

32.6 

33.3 

Wt.  % 

3.84 

2.59 

10.8 

Vol.  % 

89.5 

98,0 

95.9 

Discharge  Co- 
efficient, 

0.00312 

0.00249 

0.00292 

Cjj,  sec,"^ 

(0.01055) 

(0.01369) 

(0.00686). 

C. 

Compaction 
Density , 
gm/cc,@30  KSI 

2.60 

. 3.35 

1.53 

D. 

Burning  Rate, 

in/sec,  vs  Pg(PSIA)  q 2g 

pressed  @ 30  KSI  0.48(P/1000) 

0.21(P/1000)°'^° 

- 

E. 

Sensitivity  to 
Ignition 

Impact 

225  kg-cm 

>300  kg-cm 

42  kg-cm 

Friction 

neg  (Esso  screw) 
60“ @100  psi  (ABL) 

neg  (Esso  screw) 
75“ @300  psi  (ABL) 

neg  (Esso  screw) 
45“@100  psi  (ABL) 

Electrostatic 

Discharge 

>6J  @ 5000  V 

>6J  @ 5000  V 

0.00563  @ 5000  V 

Autoignition 

Temperature 

>400“C 

>400“C 

229*C 
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lAULK  3.12 


COMBUSTION  PRODUCTS  OF  CANDIDATE  PROPELLANTS 


Composition  3.3.4 


Specie 

Chamber 

Throat 

CT. 

0.00028 

0.00011 

H 

0.01575 

0.00597 

H2 

1.28937 

1.29426 

NH3 

0.00016 

0.00012 

N2 

0.14712 

0.14714 

OH 

0.00001 

0.00000 

ZRN8 

0.00155 

0.00000 

ZRCL4 

0.00001 

0.00000 

NK2 

0.00002 

0.00001 

ZRCL3 

0.00003 

0.00000 

HCl 

0.03363 

0.03391 

H20 

0.00328 

0.00321 

ZRCL2 

0.00001 

0.00000 

ZR028 

0.74496 

0.00004 

T "K 

2799 

2479 

$ - Condensed  Specie 


Composition  3.8.1 


Specie 

Chamber 

Throat 

CL 

0.00020 

0.00008 

HCl 

0.02487 

0.02532 

H2 

0.82210 

0.82506 

NH3 

0.00004 

0.00003 

N2 

0.01249 

0.01250 

ZRCL4 

0.00003 

0.00001 

ZR028 

0.54281 

0.00010 

H20 

0.00129 

0.00110 

WCL2- 

0.00004 

0.00002 

ZRCL2 

0.00002 

0.00000 

U8 

0.16356 

0.00003 

H 

0.00962 

0.00368 

ZRCL3 

0.00007 

0.00001 

ZRN8 

0.00050 

0.00000 

T "K 

2799 

2481 

Composition  3.8. 

8 

Specie 

Chamber 

Thro  fit 

B 

0.00011 

0.00002 

BCL2H 

0.00052 

0.00046 

BH2 

0.02901 

0.01109 

BN 

0.00003 

0.00001 

BOH 

0.00106- 

.0.00052 

B20 

0.00044 

0.00012 

CL 

0.00078 

0.00061 

HB02 

0.00074 

0.00067 

H2 

3.55937 

3.58573 

NH2 

0.00003 

0.00001 

WCL2 

0.00015 

0.00012 

ZRCL3 

0.00711 

0.00410 

ZRO 

0.00011 

0.00001 

BN$ 

1.26001 

0.04158 

ZR02$  ' 

0.35327 

0.00525 

BCL2 

0.00039 

0.00027 

BH 

0.00017 

0.00004 

BH3 

0.00249 

0.00100 

BOCL 

0.002040.0.00179 

BX03 

0.00015 

0.00012 

H 

0.04350 

0.02849 

ZRCL2 

0.00305 

0.00096 

ZRH 

0.00001 

0.00000 

BCl 

■ 0.00933 

0.00503 

BCL3 

0.00001 

0 .00001 

BH2N 

0.012‘34 

0.00543 

BO 

0.01106 

0.00570 

B202 

0.00422 

0.00243 

B3H303 

0.00001 

0.00001 

HCl 

0.09305 

0.11072 

H20 

0.00093 

0.00084 

NH3 

0.00023 

0.00014 

ZRCL 

0.00032 

0.00005 

ZRCL4 

0.00240 

0.00263 

w$ 

0.03831 

0.00003 

N2  ■ 

0.10551 

0.08825 

T “K  . 

2807 

2619 
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P = 1 000  psia  P = 75( 


MOLE  FRAC 


Figure  3.1 . Effect  of  Temperature,  Pressure 


psia 


P = 500  psia 


Concentration  on  Hy  Dissociation. 


Figure  3.2.  Burning  Rate  of  Comp.  3.3.4. 
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BURNING  I^VTE,  IPS,  X 10 


40.879%  H,WO^.  + 42.968%  ZrH^  +,  13.153%  Zr(OH),.  + 3%  NH,.C1J), 


The  propellant  composition  has  excellent  safety  properties.  Because  of  the  use 
of  NH^NO^  and  Zr(OH)^,  the  composition  mxist  be  processed  in  a low  humidity 
facility  (<50%  relative  humidity  at  70®F  is  satisfactory).  Because  of  the  phase 
change  exhibited  by  NH^NO^,  the  propellant  cannot  withstand  temperature  cycles 
of  -40/+140®F  without  loss  of  some  structural  strength. 

3.4.3  Properties  of  Composition  3.8.1 

This  composition  is  the  best  of  the  thermite  systems  (in  terms  of 
yield  and  purity)  comprised  of  off-the-shelf  ingredients.  An  excellent  clinker 
is  formed  on  combustion  at  all  pressures.  The  total  gas  yield  of  the  com- 
position is  only  2.59%  by  weight,  the  remainder  being  and  W(c)  . 

This  propellant  also  has  excellent  safety  properties.  The  use  of  Zr(OH)^ 
requires  dry  processing  and  handling  facilities  but  presents  no  major  problem. 
The  composition  can  be  temperature  cycled  with  no  adverse  affects. 

3.4.4  Properties  of  Composition  3.8.8 

This  composition  is  one  of  the  best  of  any  of  the  systems  evaluated, 
in  terms  of  yield  and  purity.  It  contains  HDB,  and  is  the  result  of  an 
extensive  laboratory-scale  search  for  ways  and  means  to  use  the  compound  in 
propellants.  It  proved  to  be  a difficult  goal  to  obtain,  because  the  material 
on  hand  was  from  various  production  lots  and  this  seemed  to  affect  propellant 
properties  to  a great  extent.  The  propellants  studied  proved  to  be  sensitive 
to  ignition  by  mechanical  means  and  difficult  to  burn  uniformly  over  a wide 
pressure  range, and  difficult  to  store  and, use,  even  in  a dry  environment. 

Initially,  work  with  the  Composition  3.9.8  (WO^/Zr/HDB/NH^ClO^) 

indicated  unstable  combustion  when  1/2  inch  diameter  pellets  pressed  at  30  KSI 

were  burned  at  1500  psi.  At  combustion  pressures  (Pc)  <1000  psi,  the  pellets 

0 58 

burn  uniformly  at  a rate  r “ 2.51  (Pc/1000)  in/sec.  At  Pc  - 1500  psi, 
these  pellets  explode  when  ignited.  When  pressed  at  20  KSI,  the  pellets  ex- 
ploded when  burned  at  400  psi.  Substitution  of  H2W0^  for  WO^  (i*e.,  Oomposition 
3.8.8)  reveals  a similar  result,  although  there  was  some  evidence  of  HDB  batch 
effects  on  this  phenomenon  since  one  of  the  batches  of  HDB  gave  uniform  burns 
at  1500  psi,  while  another  exploded.  Substitution  of  hydrazine  nitrate  (N^H^NO^) 
for  the  oxidizer  also  resulted  in  explosive  burning  at  1500  psi.  These  results 
Indicate  that  the  substitution  of  oxidizer  is  not  sufficient  to  stabilize  high 
pressure  combustion  of  HDB  propellants. 
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The  substitution  of  either  ZrH2  for  Zr  and/or  NH^Cl  (or  N2H^C1) 
for  NH^CIO^  in  propellants  containing’ HDB  and  either  WO^  or  H2^®4  gives  yield 
values  in  the  range  3-5Z  by  weight  (Comps  3.8.6  and  3.8.7  and  3.9.7)  and  as  a 
result  are  not  particularly  interestina.  Neither  compound  is  energetic  and  conse- 
auently  results  in  compositions  containing  low  values  of  HDB  content  and  therefore 
low  values  of  H2  yield.  Substitution  of  A1  for  Zr  yields  combustion  products 
which  are  significantly  contaminated  with  dissociation  products  of  Al20^ 

(which  is  not  as  thermally  stable  as  Zr02) . In  these  latter  systems,  species 
such  as  AIBO2  and  BOCl  are  also  important  combustion  products. 

The  use  of  another  chloride  as  a substitute’ for  NH^CIO^  is  possible 
in  propellants  containing  N2H^N02  as  the  oxidizer.  For  compositions. containing 
other  oxidizers,  the  favored  chlorides  are  mixtures  of  NH^CIO^  and  a fuel  (e.g. , 
Zr  or  ZrH2)  stoichiometric  to  the  major  metal  oxide  (e.g.,  Zr02>.  These  are 
high-energy  composition^  that  release  HCl  as  a major  combustion  product  and  re- 
quire large  amounts  of  coolant  (HDB, etc.)  to  maintain  the  flame  temperature  at 
the  desired  value  of  2800°K.  As  noted  abdve,  the  use  of  NH^Cl  in  place  of  NH^ 
CIO^  results  in  low  values  of  H2  yield  with  either  WO^  or  H2W0^  as  oxidizers, 
because  the  NH^Cl  is  a coolant  as  well  as  an  HCl  donor  and  the  HDB  content  is 
accordingly  reduced.  The  use  of  N2H^N0^  overcomes  this  difficulty  because  of 
the  greater  energy  of  N2H^N0^/Zr  or  ZTH2  mixtures  compared  to  those  of  WO^  or 
HzWO^/Zr  or  ZrH2.  Thus,  the  use  of  low-energy  chlorides  requires  the  use  of  a 
high-energy  oxidizer  in  order  to  maintain  acceptably  high  values  of  H2  yield. 

An  alternate  candidate  oxidizer  is  NH^NO^.  Toward  this  end,  studies  were 
directed  toward  the  system  N2H^NO^/Zr/HDB/N2H^Cl  (composition  3.5.10).  The 
commonality  of  the  ingredients  is  conducive  to  stability  and  compatibility. 

The  composition  burns  uniformly  at  pressure  up  to  1500  psi  at  .a  rate  of  r = 

1.05  (P/1000)^*^^  ips,  when  compacted  in  1/2  inch  diameter,  steel  tubes  at  30 
KSI  (compaction  density  = 1.63  gm/cc).  Unfortunately,  the  (H2  + HCl)  gas 
purity  is  (85.7%)  which  is  too  low  to  be  of  Interest.  Also,  the  composition 
is  extremely  moisture  sensitive.  At  values  of  relative  humidity  of  40%  (at 
70°F) , the  composition  physically  degrades  through  chemical  reaction.  Under 
desiccated  conditions,  no  evidence  of  chemical  reaction  is  seen.  A slight 
Improvement  in  gas  purity  can  be  attained  by  use  of  the  stoichiometric  mixtures 
Zr(OH)^  + ZrCl^  as  the  HCl  generator  (Composition  3.5.12). 
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The  best  solution  seened  to  be  through  the  use  o£  HDB  as  a mono- 
propellant coupled  with  the  oxidizer/ fuel/ chloride  as  a second  component  in 
a bipropellant  motor  configuration.  The  propellant  composition  3.8.8  can  be 
divided  into  two  propellants,  one  comprised  of  HDB,  the  other  comprised  of 
the  remainder  of  composition  3.8.8,  l.e.,  H2W0^,  Zr,  and  NH^CIO^  (called 
Mix  22).  The  physical,  ballistic,  and  safety  properties  of  these  two  propel- 
lants are  given  in  Table  3.13  and  3.14  and  Figures  3.4  and  3.5.  The  HDB 
presses  wall,  and  Ignites  and  bums  uniformly  at  all  pressures.  The  Mix  22 
bums  with  a zero  pressure  exponent,  which  is  characteristic  of  thermites.  By 
this  technique,  the  compatibility  problems  of  chemical  reaction  and  unstable 
combustion  of  the  HDB  at  high  pressures  are  obviated.  It  should  be  noted  that 
neither  composition  forma  clinkers  on  combustion,  the  HDB  because  it  is  too 
gassy,  and  Mix  22  because  it  is  too  hot  due  to  Zr02  boiling  at  4590“K.  How- 
ever, it  should  be  noted  that  the  coinbustion  products  of  Mix  22  form  clinkers 
when  mixed  with  the  combustion  products  of  HDB  in  the  bipropellant  configuration. 
This  effect  is  described  in  Section  3.5. 3. 3 below. 

3.5  Sub-Scale  Motor  Firings 

3.5.1  General 

The  objective  of  the  sub-scale  motor  firings  was  to  determine  a 
satisfactory  separation  technique  between  gases  and  solids  for  each  of  the 
candidate  propellants.  For  these  firings,  two, sizes  of  propellant  grains  were 
•valuated,  one  being  l-5/8"ln  diameter  and  the  other  being  2.5"  in  diameter.  Both 

were  cylindrical  end-burners,  inhibited  on  the  circumferential  surface  and  one 
end.  The  nominal  grain  length  was  1".  Since  the  effectiveness  of  filtration 
(or  other  separation  technique)  governs  the  mass  discharge  rate  and  therefore 
the  combustion  pressure,  it  was  the  practice  to  use  oversized  nozzles  in  first 
shots,  with  both  filter  design  and  nozzle  throat  size  changed  by  steps  to 
approach  the  desired  value  of  combustion  pressure.  It  is  recalled  that  two 
values  for  discharge  coefficient  were  given  for  each  candidate  in  Table  3.11. 

They  are  discussed  below. 
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TABLE  3.13 

PROPERTIES  OF  HDB  AND  MIX  22 


Property 

HDB 

Mix  22 

A.  Composition,  Wt.  % 

100  HDB 

19.285  H2W0^ 

54.612  Zr 

26.103  NH.CIO, 
4 4 

B.  Theoretlcala 

1.  Flame  Temperature, 
“K 

2.  (H^/HCl)  taolat 

3.  Yield  (H2+HCI) 

1878 

4590 

2.88 

Wt.  % 

16.9 

4.39 

Vol.  % 

100.0 

41.3 

4.  Discharge  coeffi- 
cient, sec”^. 

(.00237) 

(.00554) 

(0.00526) 

(0.0101) 

C.  Compaction  Density, 
gm/cc  (3  30  KSI 

0.837 

3.18 

D.  Burning  Rate,  In/ sec, 
vs  P (psla) 

0.70(P/1000)°*^® 

2.54  @ all  values  of  P 

E.  Sensitivity  to  Ignition 

1.  Impact 

41  kg-cm 

125  kg-cm 

2.  Friction 

Neg  (Esso  Strew) 
Neg  @ 90®,  400  psl 

Neg  (Esso  Screw) 
Neg  0 75®  100  psl 

3.  Electrostatic 
Discharge 

0.0013  J @ 5000  V 

.00135  0 5200 
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TABLE  3^1^ 

COMBUSTION  PRODUCTS  OF  TWO -COMPARTMENT  HDB  PROPELLANT 


HDB 


SPECIE 

CHAMBER 

BH2N 

0.00009 

NH3 

0.00005 

BN8 

3.34872 

BH2 

0.00022 

H2 

8.37201 

BH3 

0.00019 

H 

0.00078 

N2 

0.00018 

H 

o 

1877 

THROAT 

SPECIE 

0.00001 

CL 

0.00001 

HCl 

0.00048 

H2 

0.00001 

N 

8.37296 

NH3 

0.00002 

N2 

0.00010 

OH 

0.00001 

WCL 

1609 

ZRCL 

ZRCL4 

ZR02 

ZR02® 

CL2 

H20 

NH 

• NO 
02 

WCL2 
' W02 

ZRCL2 
ZRH 

H 

H2W04 

NH2 

0 

W 

WO 

W03 

ZR 

ZRCL3 

ZRO 

W* 
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MIX  22 


CHAMBER 

THROAT 

0.06385 

0.06463 

0.0995,7 

0.10410 

0.29424 

0.29971 

0.00053 

0.00038 

0.00001 

0.00001 

0.10961 

■ 0.11005 

0.01629 

0.01252 

0.00008 

0.00006 

0.00220 

0.00169 

0.00002 

0.00002 

0.00671 

0.00547 

0.54271 

0.01171 

Ov 00006 

0.00005 

0.02176 

0.0l863 

0.00010 

0.00007 

0.00224 

0.00159 

0.00026 

0.00016 

0.02599 

0.02386 

0.00114 

0.00075 

0.00140 

0.00121 

0.00008 

0.00005 

0.29482 

0.28954 

p.  00001 

0.00001 

0.00006 

0.00004 

0.00900 

0.00673 

0.00016 

0.00010 

0.00160 

0.00110 

0.00003 

0.00002 

0.00034 

0.00021 

0.00049 

0.00045 

0.04472 

0.03515 

0.04816 

0.00313 

4608 

4398 

T "K 


For  glvftn.  values  of  propellant  aasa  burning  rate  (m)  and  conib  us  tlon 
pressure  (P^) , the  noazle  discharge  coefficient  (Cj^)  Is  used  to  size  the 
nozzle  throat  area  (A^)  by  means  of  the  continuity  equation, 


m 


c D t 


(1) 


By  continuity,  m is  also  related  to  propellant  burning  surface  area^S) , 

burning  rate  (r) , and  density  (p^)  and  to  the  gas  velocity  (v^)  and 

density  (p  ) at  the  throat; 

S 


A “ Sr.p 

p 

(2) 

- A,  Pj, 

Combining  Equation  1 and  3 yields 

(3) 

S ■ ’t  Og/I'c 

(4) 

For  adiabatic  combustion  at  an  assigned  value  of  coupled  with  Isentroplc 
flow  to  sonic  velocity,  it  is  seen  that  is  thermodynamically  derived.  For 
the  motors  of  Interest,  it  is  determined  from  only  the  gaseous  combustion  products 

which  are  ejected  from  the  motor,  since  the  condensed  species  are  assumed  to 
be  retained  within  the  motor  as  a clinker  or  a filter  cake.  For  this  case 
of  uncoupled  flow  (gas  only) , the  derived  values  of  Cj^  are  the  first-reported 
numbers  in  Table  3,11,  In  the  event  that  the  condensed  species  are  ejected 
from  the  motor  because  no  clinker  is  formed  or  because  of  failure  jof  the  filter 
system,  than  a higher  value  of  is  computed.  These  values  are  re'ported  in 


parenthesis  in  Table  3.11  and  represent  an  upper  limit  on  the  theoretical  value 
for  this  parameter.  In  this  case  (coupled  flow,  i.e.,  gas  and  solids),  con- 
siderably greater  values  of  are  computed  for  a given  value  of  A^/S.  Combining 
equation  1 and  2,  and  expressing  r ■ 

1 

(P^,  coupled  \ /Cjj»  uncoupled 

P^,  uncoupled/  coupled 


where  F ■=  mass  fraction  of  gaseous  combustion  products. 


For  compositions  3.3.4,  3.8.1,  and  3.8.8,  the  P^  ratios  are  6.26,  9.55  and 
8.03,  respectively.  Thus,  for  a given  value  of  A^/S  computed  on  the  basis 
of  uncoupled  flow  (which  is  the  desired  phenomenon) , very  high  pressure  rises 
can  occur  if  the  propellant  does  not  clink  or  if  the  filter  system  falls, 
nils  is  why^  a step-vise  approach  .was  .taken,  in  slicing  the  throat  for  the  motor 
firings . 
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3.5.2  Sub-scale  Motor  Configurations 

The  inoLor  l\ar(lw/\n>  th.nL  was  const  rue t «'»<!  w;\s  t'f  a bcavv-wal  I 

case  and  closures,  wltl\  a I'ilter  housing  located  between  the  combustion  chamber 
and  nozzle  throat.  The  intent  was  to  allow  for  maximum  flexibility  in  filter 
bed  configuration  and  propellant  mass  burning  rate.  Of  the  three  candidate 
propellants,  3.8.1  clinks  best  when  burned.  For  this  composition,  a simple 
filter  bed  comprised  of  tungsten  screens  was  evaluated.  This  configuration  is 
shown  in  Figure  3.6. 

At  the  other  extreme,  composition  3.8.8  forms  no  clinker,  and  the 
solids  content  of  the  combustion  products  is  82.7%  by  weight.  For  this  com- 
position, cyclones  or  centrifugal  separators  are  best.  Attributes  of  cyclones 
include  simplicity  of  design  and  operation,  while  centrifugals  offer  a more 
compact  design  envelop.  Re- entrainment  of  particles  is  a,  phenomenon  peculiar 
to  both  techniques,  since  high  gas  velocities  are  involved.  For  the  cyclone, 
a tangential  velocity  on  the  order  of  200  ft/sec  is  required  for  separation 
of  10-y  particles  in  the  sub-scale  envelop  (3"  dia) . For  this  reason,  either 
an  impingement  separator  or  filter  bed  located  upstream  of  the  nozzle  outlet 
may  be  required  to  remove  residual  solids.  Cyclones  and  centrifugals  can  both 
be  operated  at  high  pressure.  The  principal  drawback  of  the  cyclone  is  size, 
while  that  of  centrifugals  is  the  need  for  an  external  power  source  to  drive 
the  turbine  and  the  increased  mechanical  complexity.  It  should  be  noted  that 
centrifugation  is  an  extremely  effective  technique  for  solids  removal.  In 
previous  work  performed  with  spinning  motors,  complete  removal  of  condensed 
combustion  products  (Al20^  (A) , C(c),  and  B(£)  was  obtained  for  a 6-inch  dia- 
meter motor  spinning  at  900  rpm  (centrifugal  force  -70  g's).  The  centrifugal 
effect  for  a full-scale  motor  probably  can  be  attained  by  either  spinning  the 
motor,  as  above,  or  (more  conventionally)  by  spinning  the  combustion  products 
by  means  of  a rotoclone,  or  cinder  fan.  For  compositions  3.8.8  and  3.3.4,  a 
cyclone  separation  chamber  was  employed,  primarily  because  of  design  simplicity 
and  ease  of  use  at  high  pressure.  This  motor  configuration  is  shown  in 
Figure  3.7.  The  cyclone  effect  is  imposed  on  the  combustion  products  by  the 
deflector  plate  located  above  the  propellant  burning  surface  to  direct  the 
flow  tangentially  along  the  motor  wall.  A re-entrant  nozzle  tube  was  employed 
to  minimize  the  effects  of  re-entrainment. 
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Figure  3.7.  Motor  in  Cyclone  Chamber. 


3.5.3 


Subscale  Motor  Firing  Results 


3. 5. 3.1  eomposltlon  3.8.1  (H^W0^/ZrH^/Zr(0H)^/NH^C10^) 

The  simplest  motor  configuration, that  Is  shown  In  Figure  3.6,  was 

used  for  these  firings.  The  propellant  forms  excellent  clinkers  when  burned, 

so  that  the  only  filter  arrangement  was  comprised  of  3 timgsten  screens  in 

series  located  In  the  filter  housing.  The  screens  are  1.60"  diameter  x 4 mil  x 

24  mesh,  spaced  1/4  inch  apart  in  series.  Firing  results  are  given  In  Table  3.15 

Initisl  £iLirlinss  wst*©  cv0«rs!L2wd  i®  ^ 

® “ “ ‘c  " ‘ " 

residue  weight  Increases.  Very  little  ash  was  found  on  the  screens  Indicating 
a clean  gas  situation.  Burn  rates  were  difficult  to  measure  from  the  motor 
firings  at  the  lower  values  of  D^,  because  the  ratio  of  motor  volume  to  Is 
very  large  (motor  volume  - 350  cc) , and  the  pressure  decay  times  are  significant 
fractions  of  the  burn  time.  For  no  heat  losses  the  rate  at  which  pressure  in 
the  motor  chamber  decays  to  ambient  la  derived  from  the  Perfect  Gas  Law  and  the 
continuity  equation: 


where 


dP  RT  dn 

dt  V dt 


(MW)  . 


P = Pressure  In  the  chamber 
T = Temperature  of  the  gas 
n = Number  of  moles  of  gas 


t =»  Time 
R = Gas  constant 
V = Volume  of  chamber 
MW  “Molecular  weight  of  gas 


For  firings  10-5,  -6  and  -7,  the  total  bum  times  were  6.26,  6.12  and  6.52  sec, 
respectively;  and  the  computed  decay  times  are  0.22,  0.59  and  1.44  sec,  respec- 
tively; resulting  in  computed  burn  times  of  6.26  - 0.22  = 6.04,  5.53,  and  5.08 
sec,  respectively.  The  tabulated  burning  rates,  which  are  plotted  In  Figure  3.3 
are  based  on  the  computed  times  and  are  seen  to  agree  very  well  with  strand  data 
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TABLE  3.15 
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Sub scale 

Motor  Firing  Results:  Composition  3.8.1  (H2W0^ 

/ZrH2/Zr(0H) 

^/NH^CIO^) 

Firing 

Number 

Grain  Dimensions 
Dia(in)  x L(in)  x Wt(gm) 

Nozzle 
Throat 
Diameter 
Dt,  in 

Filter 

Configuration 

Pp.  psia 

Average* 
Bum  Rate 
ips 

Residue, 

% of 

Theoretical 

Remarks 

10-1 

2.512  X 

0.586  X 122.6 

0.402 

3 W screens 

Ambient 

- 

81.5 

Ambient  burn 

10-2 

Igniter 

shot 

10-3 

1.688  X 

0.845  X 98.4 

0.128 

3 W screens 

55 

- 

93.0 

3-sec  ignition 
delay 

10-4 

1.688  X 

0.860  X 99.4 

0.089 

3 W screens 

197 

- 

100.0 

See  text 

10-5 

1.688  X 

0.842  X 98.5 

0.069 

3 W screens 

200 

0.14(0.14) 

99.7 

See  text 

10-6 

1.688  X 

0.858  X 101.5 

0.048 

3 W screens 

294 

0.155(0.15) 

97.0 

See  text 

10-7 

1.688  X 

0.845  X 98.5 

0.033 

3 W screens 

631 

0.17(0.185) 

98.5 

See  text 

*Computed  from  the  quotient  web  i computed  burn  time  (see  text)  . Numbers  in  parenthesis  are  from  strand 
data  (see  Figure  3.3). 


3. 5. 3. 2 Composition  3.3.4  (NHfN0j/ZrH^/2r(0H)^/NH^ (CO^) 

Results  of  these  firings  are  given  In  Table  3.16.  In  the  first 
firing,  9-1,  the  motor  configuration  of  Figure  3.6  was  used.  The  filter  con- 
figuration was  comprised  of  3 Tungsten  screens,  1/6"  diameter  x 4 mil  x 24  mesh 
arranged  In  series  1/4  Inch  apart.  The  firing  trace  was  very  erratic  and 
characterized  by  two  pressure  peaks.  Poat-flrlng  Inspection  Indicated  the 
screens  had  been  ruptured.  It  Is  known  that  the  composition  does  not  clink 
well  at  low  pressures,  and  apparently  the  screens  filled  and  ruptured,  causing 
the  second  pressure  peak.  The  average  pressure  Is  a poor  value,  and  the 
resultant  burning  rate  Is  therefore  also  poor. 

In  the  second  motor,  also  based  on  the  3.6  configuration,  the  filter 
housing  was  outfitted  with  three  baffles  made  of  5 mil  tungsten  sheet  and  fol- 
lowed by  a tungsten  screen.  The  spacing  between  each  baffle  was  1/4  Inch.  The 
baffles  were  1.6-lnch  dla  discs  with  a 1/4  Inch  segment  removed,  each  segment 
located  on  opposite  diameters  so  that  the  flow  through  the  baffles  was  serpentine. 
Post-flrlng  examination  revealed  that  the  baffles  were  bent  concave  In  the 
direction  of  flow,  that  the  void  space  between  baffles  was  filled  with  residue, 
and  that  the  screen  had  become  filled  and  ruptured  above  the  open  segment  of 
the  last  baffle.  The  motor  pressure  had  built  up  gradually  to  a value  of 
1500  psla  0.4  sec  after  Ignition  and  dropped  Instantaneously  to  “400  psla 
and  remained  at  approximately  this  value  for  the  remainder  of  the  firing 
(0.4  sec).  Apparently,  the  rise  to  the  pressure  spike  Is  due  to  filling  the 
baffles  and  screen,  and  the  drop  to  operating  pressure  Is  due  to  failure  of 
the  screen  after  It  had  become  completely  filled.  As  above,  the  average 
pressure  Is  not  a very  good  one,  due  to  the  spike,  and  therefore  the  motor 
and  strand  burning  rates,  which  are  shown  In  Figure  3.2,  do  not  agree  well. 

These  firings  Indicate  that  the  propellant  does  not  clink  well  enough 
to  permit  the  use  of  a simple  filter  configuration.  For  this  reason,  the 
remainder  of  the  firings  were  carried  out  with  the  motor  configuration  of 
Figure  3.7,  which  Includes  a cyclone  chamber  and  re-entrant  nozzle.  Firing 
9-3  was  carried  out  with  this  motor  and  the  same  filter  configuration  of 
9-2.  Even  though  the  P^  was  lower,  the  retention  of  residue  was  considerably 
better,  and  the  agreement  between  motor  burning  rate  and  the  strand  value 
Improved.  For  this  firing,  most  of  the  residue  was  caught  In  the  cyclone 
chamber.  None  was  found  between  the  baffles  nor  on  the  screen,  both  of  which 
were  Intact  and  not  deformed  or  ruptured  as  In  the  previous  firings.  Evidently, 
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TABLE  3.16 

Subscale  Motor  Firing  Results:  Composition  3.3.4  (NH^NO^/Zr^/Zr (0H)^/NH^C10^) 


Firing 

Number 

Grain  Dimensions 
Dia(in)  x L(in)  x Wt(gm) 

Nozzle 
Throat 
Diameter 
D^,  in 

Filter 

Configuration 

Average 
Pr,  psia 

Average* 
Burn  Rate 
ips 

Residue, 

% of 

Theoretical 

Remarks 

9-1 

2.50  X 0.947  X 194.2 

0.334 

Figure  3.6;  3 W screens 

300 

0.86  (0.33) 

55.1 

Screens  failed 

9-2 

2.51  X 0.481  X 97.7 

0.250 

Figure  3.6;  3 W baffles 
1 W screen 

400 

0.60(0.36) 

67.3 

Screens  failed 

9-3 

1.688  X 0.564  x 5-. 3 

0.099 

Figure  3.7;  3 W baffles 
1 W screen 

264 

0.40(0.32) 

93.5 

See  text 

9-4 

1.688  X 1.076  X 97.22 

0.067 

Figure  3.7;  Zr02  matt 
and  W screen 

335, 

0.43(0.35) 

100 

See  text 

9-5 

1.688  X 1.061  X 96.61 

0.055 

Figure  3.7;  Zr02  matt 
and  W screen 

627 

0.50(0.42) 

98.7 

See  text 

*Computed  from  the  quotient  web  ^ bum  time.  Numbers  in  parenthesis  are  from  strand  data  (see  Figure  3.2). 


the  cyclone  chamber  is  very  effective  in  removing  the  solids  from  the  combus- 
tion products. 

The  next  two  firings  were  also  made  with  the  re-entrant  motor  and 
with  ZrO^-fiber  filter  matt  supported  by  tungsten  screens  in  the  filter  housing 
replacing  the  baffle  screen  arrangement  of  the  previous  firings.  For  these 
firings,  retention  of  the  slag  is  essentially  quantitative.  In  all  three 
motor  firings  (9—3,  —4,  and  —5),  the  motor  burning  rates  are  significantly 
greater  than  the  strand  rates  (''-205').  The  reason  for  this  is  imexpli>ined. 

It  is  possible  that  the  grain  breaks  up  near  burnout,  resulting  in  apparently 

shorter  bum  times.  This  could  be  due  to  inhibitor  failure  Which  consisted  of 
a thinly  painted  coat  of  titanium  dioxide  filled  epoxy. 

3. 5 . 3.  3 Composition  3.8.8  (H^W0^/Zr/HDB/NH^C10_^) 

As  reported  previously,  this  composition  exhibited  combustion  insta- 
bility when  burned  at  1500  psia.  The  possibility  exists  that  pressure  spikes 
occurring  at  ignition  can  trigger  explosions  instead  of  uniform  combustion,  so 
that  no  attempts  were  made  to  bum  pellets  of  this  composition  in  motors. 
Therefore,  the  bi-propellant  two- compartment  approach  was  taken,  and  motors 
comprised  of  HDB  pellets  and  Composition  22  pellets  were  made  up.  The  web 
lengths  of  the  pellets  were  sized  to  permit  equal  burning  times,  based  on  strand 
data,  and  the  relative  number  of  pellets  was  governed  by  the  required  mass 
burning  rates  (39.97%  HDB  + 60.021%  Composition  22).  The  motors  were  comprised 
of  10  pellets  of  HDB  of  dimensions  of  0.5"  diameter  x 0.2"  L x 0.52  gm  weight, 
and  one  pellet  of  Composition  22,  0.5"  diameter  x 0.87"  L x 8.9  gm  weight.  The 
pellets  were  inhibited  on  one  end  and  the  circumferential  surface  and  functioned 
as  end  burners  in  the  motor  configuration  of  Figure  3.7  (cyclone  chamber  and  re- 
entrant nozzle) . 

Three  motors  were  fired.  Because  of  the  small  amount  of  propellant, 
residue  weights  could  not  be  measured  accurately.  Firing  results  (P^  vs  D^) 
are  given  in  Table  3.17.  Thel  first  motor  was  "fired  for  effect."  The  next 
two  were  fired  with  D^  = 0.10  inches,  the  8A-2  with  a ZrO^  fiber  matt  filter 
supported  by  a tungsten  screen  and  the  second  (8A-3)  with  two  Zr02  matt  filters 
in  series,  each  supported  by  a tungsten  screen  and  separated  by  1/4  inch.  In 
the  first  firing  (8A-2)  post  firing  observation  indicated  a light  fluffy  ash 
(BN)  throughout  the  coinbustion  chamber,  cyclone  chamber,  and  re-entrant  tube. 
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TABLE  3.17 


Subscale  Motor  Firing  Results:  HDB  + Composition  22  (H2W0^/Zr/HH^C10^) 


Firing 

Number 

Nozzle 
Throat 
Diameter 
, in 

Filter 

Configuration 

Average 
Pressure 
P^,  PSIA 

8A-1 

0.237 

Motor  3.6,  screen  only 

80 

8A-2 

0.100 

Motor  3.7,  Zr02  Matt  + W Screen 

575 

8A-3 

0.100 

Motor  3.7,  2 Zr02  Matts  + W screens 

560 

The  combustion  chamber  also  contained  an  acom-sized  metallic  clinker  weighing 
3.7  gm,  which  was  the  residue  from  Composition  22.  The  filter  and  screen  were 
intact. 

The  residue  from  8A-3  was  very  similar  to  the  above,  but  the  filters 
had  been  bent.  Apparently  the  use  of  two  filters  gives  rise  to  a value  of  AP 
across  the  filters  that  is  too  high.  Since  there  was  no  difference  in  the 
amount  of  residue  on  the  filters  between  the  second  and  third  test,  one  filter 
would  prove  sufficient. 

The  motor  firing  traces  for  firing  8A-2  and  8A-3  are  shown  in  Figure 
3.8a  and  3.8b.  As  can  be  seen  from  the  traces,  the  chamber  pressure  is  essentially 
constant.  Thus,  even  though  the  HDB  la  tempermental  when  combined  with  other 
ingredients,  in  a two-compartment  motor  the  ballistic  properties  of  this  system 
are  excellent. 

Theoretical  calculations  of  the  propellant  combustion  products  have 
been  made  and  have  been  compared  with  experimental  values.  The  experimental 
procedure  used  in  the  analysis  was  to  combust  small  pressed  propellant  samples 
in  a stainless  steel  bomb.  The  Interior  of  the  bomb  and  all  exposed  metal 
surfaces  were  sandblasted,  cleaned  with  trichloroethylene,  then  teflon  coated. 

The  teflon  coating  was  primarily  used  to  prevent  hbsorption  of  HCl  on  cold 
metal  surfaces  - a phenomenon  observed  under  previous  programs.  Before  and 
after  teflon  coating,  the  bomb  was  baked  out  in  an  oven  overnight  at  800°F 
to  remove  any  absorbed  gases.  All  fittings,  valves  and  transfer  tubing  was 
constructed  of  teflon  or  a similar  HCl  unreactive  plastic. 

The  bomb  was  purged  a minimum  of  5 times  with  argon  and  vacuum 
pumped  between  purges  to  ensure  no  contamination  with  residual  air.  Bomb 
volume  was  exactly  measured  via  displacement  and  the  propellant  weight  recorded. 
Ignition  was  carried  out  by  a hot  wire  technique,  to  avoid  any  introduction 
of  contaminating  gases,  under  one  atmosphere  of  argon.  Upon  cooling,  a plastic 
syringe  sampled  gases  through  a septum  fitted  directly  to  the  bomb  and  trans- 
ferred the  gases  directly  to  a calibrated  gas  chromatograph  for  analysis. 

Transit  time  and  residence  time  in  the  bomb  were  minimized.  The  results  of  a 
gas  chromatograph  of  the  sample  gave  results  that  essentially  agreed  with  pre- 
dicted thermochemistry  calculations. 
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PRESSURE  (psia) 
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^urn  Time  = 0.37  sec  versus  0.34  from  strands 
P = 575  psia 

Grain  Configuration:  10  HDB  pellets  0.2  in.  L X 0.495  dia. 

1 Comp.  22  pellet  0.87  in.  L X 0.495  dia. 

All  pellets  end  burners  burning  on  one  end 
Total  propagation  weight  = 14.1  gm 

Motor  Configuration:  Re-entrant  nozzle  with  cyclone  chamber 

One  Zr02  fiber  mat  filter  with  tungsten  screen 
Dj  = 0.1  in. 

Post-firing  Configuration:  Comp.  22  clinker  in  combustion  chamber.  Bn  is  uniformaly 


distributed  in  combustion  cyclone,  and  re-entrant  chambers. 
Filter,  filter  cake  fBnl  and  screen  are  intact. 


TIME  (sec) 

Figure  3 . 8b 


FIGURE  3 . 8 


3.6 


SUMMARY 


The  objective  of  the  solid  propellant  technology  effort  was  to 
select  and  characterize  three  solid  propellant  candidates.  Selection  was 

f 

based  on  the  results  of  calculations  of  theoretical  flame  temperature  and 
combustion  product  composition.  For  these  calculations,  propellant  Ingre- 
dient selection  was  constrained  by  state-of-the-art  solid  propellant 
chemistry.  Successful  sub-scale  motor  firings  (1/4  lb  nominal  propellant 
weight)  were  carried  out  for  all  three  candidates.  Design  factors  for  fil- 
tration and  particle  separation  techniques  for  each  candidate  were  delineated. 

Of  the  three  candidates,  the  first  represents  the  best  clinker  pro- 
pellant. This  propellant  Is  comprised  of  H^WO^,  ZrH2,  Zr(OH)^,  and  NH^ClO^, 
has  a flame  temperature  of  2799°K  and  an  (H2  + HCl)  yield  of  2.59%  by  weight 
at  98.0%  gas  purity  and  H2/HC1  ■>  32.6  (molar).  A simple  filter  comprised 
of  a tungsten  screen  was  found  to  be  satisfactory  for  retaining  entrained 
solids  In  the  combustion  products.  Typical  clinker  weight  >99%  of  theoreti- 
cal values  were  measured  for  this  propellant.  The  second  propellant  repre- 
sents the  best  attainable  from  state-of-the-art  ingredients.  The  propellant 
Is  comprised  of  NH^NO^,  Z*'H2j  Zr(OH)^,  and  NH^CIO^,  has  a flame  temperature 
of  2801“K,  and  an  (H2  + HCl)  yield  of  3.84%  by  weight  at  89.5%  gas  purity 
and  H2/HCI  = 38.2  (molar).  Although  this  propellant  does  not  form  clinkers 
as  well  as  the  H2W0^  system,  clinker  weights  >99%  of  theoretical  were  attained 
through  the  use  of  a cyclone  chamber  and  a filter  bed  comprised  of  Zr02  fiber 
matt.  A similar  particle  separator/filtration  technique  was  used  satisfactorily 
for  the  third  propellant  which  Is  the  best  attainable  from  use  of  advanced 
propellant  ingredients.  This  propellant  Is  comprised  of  H2W0^,  Zr,  HDB  and 
NH^CIO^,  has  a flame  temperature  of  2808®K,  and  an  (H2  + HCl)  yield  = 10.8% 
by  weight  at  95.9%  gas  purity  and  H2/HCI  =■  33.3  (molar).  In  this  propellant, 

HDB  (hydrazine  diborane)  is  the  codlant,  HDB  is  a compound  that  is  easily 
synthesized  and  is  available  on  special  order.  The  motor  firings  for  this 
propellant  are  especially  unique  because  a bipropellant  grain  configuration 
was  used  in  the  motors  due  to  the  reactivity  during  combustion  of  HDB  with 
NH^CIO^.  One  grain  was  comprised  of  H2W0^,  Zr,  and  NH^CIO^  while  the  other 
was  comprised  of  only  HDB.  The  firings  yield  uniform  combustion  and  pressure 
time  traces,  and  the  filtration  system  withstood  the  firing  conditions  with- 
out damage. 
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The  high  values  of  yiilfl  and  purity  of  (H2  + HCl)  offered  by  use 
of  HDB  in  propellants  leads  to  the  conclusion  that  additional  work  with 
this  and  possibly  other  HDB  propellants  is  warranted.  Synthesis  of  additional 
HDB  is  needed  to  permit  propellant  evaluation  with  a single  lot  of  HDB,  elim- 
inating batch  effects.  Scale-up  firings  should  be  made,  with  the  emphasis  on 
BN(c)  separation.  The  combustion  products  of  H2W0^/Zr/NH^C10^  form  a good 
clinker  when  admixed  with  the  autefoombustlon  products  (H2,  BN)  of  HDB,  and  for 
this  reason  the  use  of  a cyclone  or  centrifugal  separator  for  only  HDB  combus- 
tion products  minimizes  the  envelope  and  mechanical  Impact  on  the  system  as 
well  as  the  heat  sink  effect  of  the  filter  hardware. 
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4.0 


PROOF  EXPERIMENT  AND  LASER  SYSTEM  CHARACTERISTICS 


The  end  objective  of  the  work  described  in  this  report  is  to  develop 
an  H2-HCI  laser  device.  Prior  to  design  of  such  a device,  it  Is  anticipated 
that  a laser  proof  experiment  will  be  performed.  The  objective  of  the  proof 
experiment  Is  to  experimentally  determine  laser  performance  at  conditions  that 
are  of  practical  Interest.  The  current  phase  is  concerned  with  determining 
a low  risk  approach  for  the  proof  experiment  and  to  develop  the  techniques 
that  will  be  used  in  design  of  the  proof  experiment  hardware. 

An  early  objective  of  this  task  was  directed  at  obtaining  perfor- 
mance maps  as  a function  of  gas  composition,  chamber  conditions,  and  nozzle 
contour.  The  one-dimensional  harmonic  oscillator  program  was  used  to  generate 
these  maps  using  the  vibrational  kinetics  for  the  H2/HCI  system  listed  in 
Table  2.2.  These  calculations  provided  (i)  a limited  performance  optimiza- 
tion, (11)  a base  line  performance  for  conditions  expected  to  be  achievable 
by  propellant  technology,  and  (ill)  some  effects  on  H-atom  performance  de- 
gradation. 

Some  preliminary  performance  sensitivity  studies  on  stagnation  and 
gas  expansion  conditions  has  been  described  in  Phase  I.  It  was  shown  that 
best  performance  was  obtained  at  highest  stagnation  temperature  (neglecting 
atom  effects).  However,  it  was  also  determined  that  there  was  an  optlmtnn 
value  of  p^  h^^  where  p^  in  stagnation  pressure  and  h is  the  nozzle  throat 
height.  This  effect  is  similar  to  that  observed  in  the  CO2/N2  GDL  and  is 
easily  explained  in  terms  of  the  H2  vibrational  energy  freezing  criterion. 

In  addition,  it  was  shown  that  maximum  performance  occurred  at  lowest  HCl 
concentrations . This  latter  conclusion  may  be  expected  to  be  modified  in 
favor  of  larger  HCl  concentrations  in  the  A.H.O.  model  and  when  cavity 
power  extraction  is  included  in  the  modeling. 

Some  additional  calculations  illustrating  the  points  mentioned 
above  are  given  in  Figures  4.1  and  4.2.  It  is  seen  that  maximum  available 
energy  occurs  at  [HCl]  'i<2%  for  the  smaller . throat  height  and  * 50  atm. 

At  the  higher  stagnation  pressure  the  throat  height  has  a large  effect  on 
the  performance  whereas  at  the  lower  chamtber  pressure  the  performance  is 
less  sensitive.  This  effect  is  similar  to  that  which  occurs  in  the  CO2 
GDL  where  performance  drops  off  at  large  valxjes  of  the  parameter  p^  ^th' 
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AVAILABLE  POWER  (kj/lbm) 


SMALL  SIGNAL  GAIN  (percent/cm) 


AVAILABLE  POWER  (kj/lbm) 


Figure  A. 2 
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SMALL  SIGNAL  GAIN  (percent/cm) 


Some  calculations  were  also  performed  using  various  assumed  nozzle 
contours  (or  gas  density  distributions).  The  assumed  density  distributions 
are  shown  In  Figure  4.3.  Results  of  this  limited  study  Is  Table  4.1.  it 
Is  clear  that  some  Improvement  In  performance  Is  possible  by  suitably  tailor- 
ing the  nozzle  contour  to  maximize  H2  -*■  HCl  V->V  pumping.  This  study  was  not 
pursued  further  pending  resolution  of  the  basic  kinetics  of  the  H2/HCI  system. 

Another  technique  was  used  to  estimate  the  optimum  area  ratio  for 
GDL  performance.  A calculation  was  performed  comparing  the  characteristic 
reaction  times  for  the  various  vibrational  energy  transfer  processes.  An 
Isentroplc  flow  was  used  to  calculate  the  thermodynamic  state  associated  with 
each  expansion  ratio  starting  from  some  assumed  stagnation  conditions.  The 
characteristic  time  of  the  reaction  Is  equivalent  to  the  vibrational  relaxa- 
tion time  for  that  process  at  the  particular  temperature  and  density  that 
result  from  the  expansion.  An  example  of  such  a comparison  Is  shown  In  Figure 
4.4.  Also  In  Figure  4.4.  Is  a dashed  line  that  Indicates  the  time  for  an 
expanded  gas  sample  to  traverse  a one  meter  long  cavity.  For  an  efficient 
device  we  wish  the  V->V  transfer  reaction  time  to  be  short  compared  to  the 
cavity  flow  time  and  In  turn  shorter  than  all  V->-T  deactivation  processes. 

It  can  be  seen  that  In  Figure  4.4  this  criterion  occurs  generally  for  A/A* 

< 30.  Calculations  like  this  were  performed  for  a variety  of  concentration 
ratios  and  stagnation  conditions  to  develop  Insight  Into  optimum  expansion 
ratio  for  H2/HCI  vibrational  pumping. 

The  results  of  this  analysis  Is  shown  In  Figure  4.5.  The  available 
energy  In  the  HCl  at  a distance  of  100  cm  from  the  throat  Is  given  as  a func- 
tion of  area  ratio  for  several  HCl  concentrations.  Energies  greater  than 
100  kJ/lb  are  obtained  for  HCl  mole  fractions  of  0.01  - 0.02  for  A/A*  50. 
This  result  Is  comparable  to  the  previous  results  In  the  sense  that  It  gives 
a rough  estimate  for  the  nozzle  expansion  ratio.  It  Is  also  apparent  from 
this  limited  optimization  study  that  the  stagnation  temperature  would  be 
lowered  to  2800°K  where  the  propellant  technology  Is  more  practical  with  no 
performance  degradation.  It  was  recognized  In  the  Phase  I feasibility  study 
that  H-atoms  which  occur  at  the  high  stagnation  temperature  and  freeze  during 
the  rapid  expansion  can  seriously  degrade  calculated  performance  by  rapidly 
deactivating  HCl.  During  this  program  a series  of  calculations  was  performed 
at  successively  lower  stagnation  temperature.  Of  course,  It  was  expected 
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TABLE  4.1 


Initial  conditions : T 

Contour 

old 

I 

II 


=>  3000 “K,  P atm,  [HCl]  = 0.02 
o 


Available  Energy 
kJ/lb 


93 

109 

116 
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CHARACTERISTIC  TIME,  t.  sec 


Figure  4.4  Characteristic  Times  for  Various  Energy 

Transfer  Processes  as  a Function  of  Expansion  Ratio. 
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AVAILABLE  ENERGY,  Kj/lb. 


Figure  4.5  Optimized  performance  map  for 
baseline  condition  of  Tg  = 2800°K 
®-iid  Pg  = 40  atm. 
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that  specific  energy  would  decrease,  but,  at  the  lower  T , H-atoma  would  also 

o 

decrease  and  have  less  effect  of  performance. 

The  method  employed  was  to  systematically  modify  the  nozzle  contour 
and  observe  the  change  in  the  vibrational  temperature  of  HCl  versus  axial 
distance.  At  temperatures  of  2400°K  and  2300“K,  the  specific  power  computed 
is  45  kJ/lb  and  25  kJ/lb,  whereas  the  gain  is  .55%/cm  and  .35%/cm,  respectively. 
Although  these  values  are  not  as  attractive  as  the  values  obtained  at  the 
higher  temperatures,  the  effect  of  H-atoms  at  2300“K  is  estimated  to  be  neg- 
ligible. Thus,  from  the  point  of  view  of  a low  risk  proof  experiment  a stag- 
nation temperature  of  2300“K  appears  to  be  a good  choice. 

These  estimates  use  the  equilibrium  throat  value  of  H-atoms.  The 
area  ratio  of  the  resulting  nozzle  giving  the  best  performance  was  15.  For 
a throat  height  of  .01  cm,  the  overall  length  of  the  nozzle  is  2 cm.  Since 
the  nozzle  has  such  a low  area  ratio,  a slight  Increase  in  gain  can  be  ob- 
tained without  significant  degradation  of  pressure  recovery  by  placing  an 
additional  expansion  downstream  of  the  nozzle.  The  vibrational  temperature 
of  HCl  for  this  case  is  shown  in  Figure  4.6  along  with  the  results  of  other 
calculations  at  different  expansion  ratios  and  chamber  pressures.  Another 
implication  of  Figure  4.6  is  that  the  performance  is  not  as  critical  on  the 
H2  vibrational  temperature  as  originally  believed.  A stagnation  temperature 
of  2300 “K  for  the  H2-HCI  mixture  gives  a frozen  vibrational  temperature  of 
2100®K  for  the  H2.  This  relatively  low  H2  vibrational  temperature  pumps  the 
HCl  to  a vibrational  temperature  of  6500“K. 

The  harmonic  oscillator  program  was  then  modified  to  include  the 
H-atom  recombination  kinetics  to  see  if  the  equilibrium  H-atom  throat  value 
could  be  obtained  xjslng  a quick  freeze  nozzle.  The  original  contour  had  a 
subsonic  approach  section  of  .1  cm  long.  This  resulted  in  essentially  no 
H-atom  recombination  in  the  nozzle.  The  subsonic  approach  section  was  in- 
creased to  1.6  cm  and  3.2  cm  as  shown  in  Figure  4.7.  The  resulting  H-atom 
concentration  versus  x is  shown  in  Figure  4.8.  The  theoretical  results 
Indicate  that  the  equilibrium  throat  value  of  H-atom  concentration  is  not 
achieved.  A similar  calculation  was  performed  for  chamber  conditions  of 
80  atmospheres  and  2800°K.  The  H-atom  concentration  is  shown  in  Figure  4.9. 
Again,  the  H-atoms  do  not  reach  the  throat  equilibrium  value  although  at  the 
large  subsonic  lengths  the  differences  are  not  large. 
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Figure  4.6  HCL  Vibrational  Temperature. 
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FIGURE  4.7.,  NOZZLE  SUBSONIC  SECTION  DENSITY  DISTRIBUTION 
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FIGURE  4.8.  H-ATOM  CONCENTRATION  IN  SUBSONIC  NOZZLE  REGION 
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FIGURE  4.9 


H-ATOM  CONCENTRATION 


4.1 


Method  of  Characteristics  Analysis 


The  harmonic  oscillator  or  anharmonlc  oscillator  analysis  Is  not 
capable  of  generating  a smooth  nozzle  contour  that  will  produce  a nozzle 
exit  flow  that  is  uniform  and  shock  free.  For  this  reason,  a method  of 
characteristics  computer  program  that  takes  into  account  the  nonequilibrium 
nature  of  the  flow  was  developed.  This  program  is  a combination  of  a method 
of  characteristics  program  and  the  harmonic  oscillator  analysis  program. 

The  modifications  required  to  the  method  of  characteristics  program  were  to 
use  the  frozen  sound  speed  defined  as 


where  p “ pressure  R “ gas  constant 

p - density 

where  C are  comprised  of  only  the  rotational  and  translational  contributions 
to  the  specific  heat.  Also,  the  velocity  u and  temperature  T that  appear  in 
the  compatibility  equations  are  determined  from  the  solution  of  the  energy  and 
momentum  equations  which  are : 


uG  dT  u dh  R E 
d^- 

The  terms  involving  Ey  represent  the  effect  of  vibrational  nonequilibrium  and 
are  the  sum  of  the  compounds  from  each  specie.  The  quantity  h is  the  static 
enthalpy  as 

h = E . + C T 
1 vi  p 

where  is  the  mole  fraction  of  specie  i. 
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The  compatablllty  equations  for  the  characteristics  solution  also  contains  a 
heat  release  term  which  accounts  for  nonlsentroplc  vibrational  energy 
changes.  The  compatibility  equations  are 


A dp  + de  + B 


. S Xi  Ci± 

2 

u^CpT 


where 


A 


B 


uy 


3 ^vl 
3t 


The  Increments  dn  and  d^  are  shown  In  Figure  4.10  and..are  the  distances  along  the 
uprunnlng  and  downrunning  characteristics,  respectively. 


The  computer  program  has  been  used  to  determine  a preliminary 
contour  for  the  Mod  VIII  nozzle  contour.  The  area  distribution  from  the 
harmonic  oscillator  program  was  used  to  determine  an  Initial  estimate  for 
the  contour.  This  was  then  used  to  determine  the  nozzle  slope  as  a function 
of  X.  This  Initial  flow  angle  distribution  Is  shovm  In  Figure  4.11.  This 
Initial  contour  produced  the  wave  diagram  In  Figure  4.12  and  caused  the  pro- 
gram to  terminate  at  about  x = .125  cm  on  the  centerline.  This  was  due  to 
the  formation  of  a shock.  The  contour  was  then  subsequently  modified  to  that 
designated  as  A In  Figure  4.11.  The  wave  diagram  for  this  case  Is  shown  In 
Figure  4.13.  As  can  be  seen  from  the  figure,  a weaker. shock  Is  now  forming 
due  to  too  rapid  turning  of  the  nozzle  wall  back  to  zero.  The  shock  Is 
weaker  than  the  Initial  contour  but  Is  still  present.  Further  reduction  In 
the  nozzle  wall  angle  would  eventually  result  In  completely  eliminating  the 
shock. 
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Figure  4.11 


Mod  VIII  Nozzle  Contour,  0 Versus 


4.2 


Botmdary  Layer  Analysis 


The  invlsctd  nozzle  contour  as  determined  from  the  method  of  charac- 
teristics analysis  has  to  be  corrected  for  viscous  effects  due  to  boundary 
layer  growth.  To  accomplish  thla»  a bovindary  layer  integral  technique  was 
used.  Due  to  the  extremely  large  flow  gradients  in  the  throat  region  of  a 
quick  freeze  gas  dynamic  laser  nozzle,  an  integral  approach  requires  considerably 
less  computer  time  than  a finite  difference  technique.  The  integral  technique 
employs  the  momentum  Integral  equation  which  is 


ii  + 9 d_  ,(P  . 6*  dUp 

dx  „2  dx  ® ® U dx 

p U e 

e e 


T. 


e e 


where  = boundary  layer  edge  velocity 

6*  = displacement  thickness 
6 = momentum  thickness 
6 = boundairy  layer  thickness 
= boundary  layer  edge  density 
T = wall  shear  stress 


There  are  three  unknowns  in  this  equation,  0,  6*  and  Tu  . The  quantities  6* 
and  0 are  defined  as 


o 


The  evaluation  of  6*  and  0 require  specifying  the  velocity  profile 
and  the  density  profile.  The  velocity  profile  selected  is  a power  law  profile 
of  the  form 


(n) 


l/N 
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The  quantity  N is  determined  using  the  empirical  data  of  Reference  40.  A curve 
fit  of  N is  given  in  Figure  4.14.  As  can  be  seen  in  the  figure,  the  value  of 
6 and  subsequently  R^0  determines  the  value  of  N which  in  turn  determines  th<=i 
velocity  profile.  This  approach  is  simpler  than  resorting  to  the  use  of  an 
additional  equation  to  determine  N and  since  the  real  test  of  the  method  must 
rely  on  experimental  data,  for  the  present  this  technique  is  probably  sufficient. 

The  density  profile  and  ^ear  stress  were  determined  using  the  results 
of  Reference  41.  The  work  of  Reference  41  reviewed  the  literature  on  compressible 
turbulent  boundary  layers.  One  of  the  results  of  Reference  41  was  to  develop  an 
expression  for  the  skin  friction  coefficient  from  which  the  wall  shear  stress 
T(j  can  be  determined.  The  skin  friction  relation  developed  has  the  form 

F R 
f F 

c 

where 


The  density  distribution  used  in  the  above  relation  is 


fi-  = 


-1 


JL+  1 + i r (y  - 1)  ^ u_  r (y  - 1)  M?  /u_\^ 

T-  2 T.  u 2 J 


where  r is  the  recovery  factor  and  has  a value  of  about  .89.  A plot  of 

versus  F „ R .taken  from  reference  41  shown  in  Figure  4.15.  These  values  are 
r6  e6 

also  tabulated  in  the  reference  and  are  given  in  Table  4.2.  The  two  sets  of  data, 
however,  do  not  agree.  In  order  to  eliminate  this  problem  a test  case  for 
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Figure  4JL4  Variation  of  N with  Reg  for  Axisymmetric  Air  and 
Nitrogen  Nozzle  Data. 
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Figure  4.15  Comparison  Between  Theoretical  and  Experimental  F^Cf  Versus 

0,  Experiments,  Adiabatic-,  O,  Experiments  with  Heat  Transfer;  Theory. 


TABLE  4.2 


BOUNDARY  LAYER  CHARACTERISTICS  * 


F c. 
c f 

^Rd^e0 

F c, 
c f 

R(j  e0 

0*0010 

2*878  X 10^ 

0*0060 

233*0 

0*0015 

3*955  X 10^ 

0*0065 

177*6 

0*0020 

5*425  X 10^ 

0*0070 

140*4 

0*0025 

1*386  X 10^ 

0*0075 

114*4 

0*0030 

5030 

0*0080 

95*62 

0*0035 

2283 

0*0085 

92*49 

0*0040 

1208 

0*0090 

70.91 

0*0045 

716*0 

0*0095 

62*55 

0*0050 

462*3 

0*0100 

55*87 

0*0055 

319*4 

0*0105 

50*46 

*Reference  41  - "The  Drag  of  a Compressible  Turbulent  Boundary 
Layer  on  a Smooth  Flat  Plat  with  and  without  Heat  Transfer," 

D.B.  Spalding  and  S.W.  Chi,  Journal  of  Fluid  Mechanics,  Vol.  18, 
Part  I,  January  1964. 
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air  was  run  that  employed  the  analyils  of  reference  42.  The  results  of  this 
case  ark  shown  in  Figure  4.16..  ^e  curve  fit  for  employed  in  the  present 
technique  that  gave  the  best  agreement  with  the  method  of  reference  42  was 


.0462 


(^rd“ee) 


-.283826 


/ F. 


The  technique  was  then  used  to  determine  the  boundaiy  layer  properties 
for  the  Mod  VIII  nozzle  contour.  The  results  of  this  calculation  are  shown  in 
Figure  4.17.  The  results  Indicate  that  the  boundary  layer  displacement  thick- 
ness is  relatively  small  and  represents  only  a small  percentage  of  the  nozzle 
exit  flow  area.  It  must  be  remembered  that  these  results  are  valid  only  within 
the  limits  of  the  analytical  model  employed  but  represent  .a  good  first  estimate 
of  the  viscous  effects. 
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Figure  4.16 
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Figure  4.17  Nozzle  Boundary  Layer  Properties. 
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4.3 


Propellant  Impurity  Effects 


During  the  solid  propellant  gas  generator  development,  it  was  anti- 
cipated that  various  species  other  than  H2/HCI  and  their  atoms  might  also  be 
produced  in  small  concentrations.  A list  of  potential  contaminant  species  is 
given  in  Table  4.3.  The  effect  of  these  Impurity  species  on  predicted  laser 
performance  is  important,  since  it  could  well  determine  the  technical  direction 
of  any  further  propellant  development  work  undertaken. 

Most  of  the  impurity  species  listed  in  Table  4.3  are  metallic  halides, 
oxides , or  hydroxides  with  relatively  high  boiling  points . Because  they  are 
not  easy  to  produce  in  the  gas  phase,  there  is  no  experimental  data  on  their 
efficiency  of  vibrational  relaxation  on  HCl  or  H2*  The  fundamental  vibrational 
frequencies  of  these  refractory  species  are  quite  small  compared  to  those  for 
and  HCl.  Therefore,  it  is  anticipated  that  these  impurities  will  not 
act  as  V ->  V deactivators . On  the  other  hand,  the  efficiency  of  a V T 
deactivation  might  reasonably  be  expected  to  vary  between  a probability  per 
collision,  P,  of  10  ^ to  10  For  Illustration,  we  will  assume  a ’i^alue  of 

_3 

10  which  is  also  similar  to  the  probability  for  the  HCl/HCl  V->-T  process  at 
T = 300°K.  Therefore,  for  impurity  specie  concentrations  <<[HC1],  the  effect 
of  the  contaminant  should  be  insignificant.  To  allow  for  the  possibility  of 
a larger  probability  per  collision,  a more  conservative  estimate  of  [impurity] 
<0.1  mole  percent  is  suggested.  This  appears  to  be  generally  within  the  capa- 
bility of  the  current  propellant  technology. 

Of  the  species  listed  in  Table  4. 3 only  two  possess  fundamental  fre- 
quencies close  to  those  of  HCl,  e.g.,  H20,and  N2.  Some  data  do  exist  on  the 

vibrational  coupling  of  N2  to  HCl,  so  it  was  decided  to  do  a more  quantitative 
estimate  of  the  effect  of  these  species  on  laser  performance. 

H„0  is  a special  case.  It  is  known  that  H_0  deactivates  HCl  very 

35  ^ 

rapidly  , although  no  temperature  dependence  is  available.  If  the  large 

cross-section  for  deactivation  persists  to  elevated  temperature,  the  effect  of 

H2O  on  laser  performance  would  be  veiry  similar  to  that  of  H-atoms.  For  this 

reason,  the  concentration  of  H2O  was  kept  to  a low  value  in  the  propellant 

formulations  selected. 
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TABLE  4.3 


POTENTIAL  GAS  CONTAMINANTS 


SPECIES 

BF 

BF2 

BF3 

ZrCl 

ZrCl2 

ZrCl3 

ZrCl^ 

AlCl 

AICI2 

AICI3 

AlOCl 


AI2O 

AlOH 


H2O 

N. 


VIBRATIONAL  FREQUENCIES 
(cm 

1410 

1120,  470,  1450 
881,  670,  1463,  480 
437 

346,  92,  461 
490,  185,  333,  131 
376,  95,  418,  110 
480 

430,  210,  540 
371,  146,  185,  610 
450,  350,  900 
715,  238,  994 
1000,  1300,  3600 
3652.  1595,  3756 
2331 
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vibrational  rate  constants  have  been  assembled  for  N„/H„/HC1  molecu- 
lar  systems  from  the  available  data  * ’ and  are  listed  in  Table  4.4.  Figure 

4.18  shows  the  result  of  a calculation  of  available  energy  as  a function  of  N2 
mole  fraction  for  the  assumed  base  line  condition.  The  results  of  this  calcu- 
lation indicate  that  2-3%  of  N2  is  acceptable  without  appreciably  affecting 
performance.  A major  part  of  the  performance  change  is  due  to  the  Increased 
effective  mass  of  gas  as  N2  is  added.  Some  of  the  loss  can  be  recovered  through 
improved  gas  yields  made  possible  by  permitting  N2  formation.  Thus,  the  gas 
generators  selected  are  sufficient  for  good  performance;  however,  a more  de- 
tailed assessment  of  N2  must  await  better  kinetic  data  or  shock  tunnel  gain 
experiments . 


4.4  Power  Extraction 

The  initial  feasibility  calculations  were  all  based  on  the  concept 

of  "available  energy,"  l.e. , the  amount  of  potentially  abstractable  energy 

existing  in  HCl  at  some  selected  position  in  the  cavity.  To  provide  a more 

realistic  analytic  capability,  power  extraction  models  based  on  the  "gain 

equals  loss"  approximation  were  developed  and  became  operational  for  both  the 

H.O.  and  A.H.O.  codes.  These  analyses  for  power  extraction  are  based  on  the 
35 

Rlgrod  approach  employing  the  concept  that  lasing  can  occur  when  the  gain 
exceeds  the  single  pass  cavity  loss.  The  photon  amplification  condition  in 
the  laser  cavity  can  be  expressed  as: 

^(2g^  ■ ^ Jl  ) > [R  (1  - (32) 

where  g^  ^ is  the  gain  for  vibrational  transition  v^v  - 1,  Ji  is  the  cavity 
length  between  mirrors  of  "effective"  reflectivity  R (which  includes  all  use- 
less losses),  and  is  the  cavity  output  coupling.  Equation  (32)  can  be 
rewritten  as  the  lasing  condition 


In  [R  (1  - L )] 
c 

2 a 


(33) 


The  gain  is  determined  from  the  level  populations  which  are  in  turn  determined 
by  the  kinetic  equations . 

The  effect  of  stimulated  emission  is  Incorporated  into  the  differen- 
tial equations  for  the  level  populations  by  adding  to  the  kinetic  terms  the 
effects  of  lasing  to  and  from  the  particular  level;  namely 
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TABLE  4.4 


VIBRATIONAL  KINETICS  FOR  N^/H|,/HCl  SYSTEMS 


3 -1  -1 

(Rate  constants  in  units  of  cm  - particles  sec  ) 


V T Process  (26) 

HCl*  + N2  HCl  + N2 

kf  = 5.53  X 10  T exp  (-109.62  T cm^-molecule  ^-sec  ^ 

H2*  + N2  ^ H2  + N2  (27) 

kf  - 5.44  X 10“^°  exp  (-114.91  t“^^^) 

N2*  + HCl  -»■  N2  + HCl  (28) 

(not  known,  but  unimportant  when  H2  HCl) 

N2*  + H2  N2  + H2  (29) 

k^  = 1.78  X 10~^  exp  (-172.69  T"^^^) 

N*  + N2  N2  + N2  (30) 

-7  -I/3 

k^  - 8.53  X 10  exp  (-273.10  T 

V ->■  V Process 

HCl*  + N2  -»■  HCl  + N*  (31) 

k^  = 8.10  X lO"^^  T 

(The  H2/N2  V V process  is  considered  as  unimportant  and  thus  neglected) 
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Figure  4.18.  Effect  of  N£  on  Performance. 
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kinetic 


(34) 


d n d n 

V ^ V 

dt  dt 


j.  V „v  V - 1 

®v  + 1 + 1 " ®v 


Is  the  cavity  flux  determined  from  the  Imposition  of  Equation  (33), 
This  model  assumes  rotational  equlllbrlvun.  For  a given  vibrational 
level,  only  the  rotational  line  of  largest  gain  Is  allowed  to  laser,  but 
simultaneous  lasing  on  different  vibrational  levels  Is  allowed.  In  the  GDL 
case,  the  total  power  extracted  la  calculated  from  the  expression: 


E - £ J’  hv;;  ■ ^ ^ A/A*  4r  ^ »5) 

V p u 

* * V - 1 

where  p and  u are  the  gas  density  and  velocity  at  the  throat,  and  hv^ 

Is  the  energy  of  the  transition.  The  results  of  a calculation  using  this 

model  are  Illustrated  In  Figures  4.19  - 4.20. 

Figure  4.19  shows  the  calculated  small  signal  gain,  g^ , determined 
for  the  given  conditions  from  the  A.H.O.  kinetics  model.  The  maximum  gain 
occurring  on  vibrational  levels  v = 1 - 7 Is  shown  as  a function  of  distance 
downstream  of  the  nozzle  throat.  Assuming  an  optical  cavity  condition  of 
R = 0.98  and  = 0.2  a critical  gain  can  be  calculated  from  Equation  (33) 
and  Is  shown  on  Figure  4.19  as  the  dashed  line.  Lasing  should  occur  for  x ^ 
10  cm  downstream  of  the  throat. 

Figure  4.20  shows  the  same  calculation  as  Figure  4.19  Including 
power  extraction  at  distances  greater  than  15  cm  downstream.  The  effect  of 
lasing  Is  seen  to  surpress  the  gain  to  values  around  the  critical  value. 

Only  a few  gain  histories  are  plotted  on  Figure  4.20  for  clarity.  The  ob- 
served oscillatory  behavior  Is  due  to  rotational  line  shifting  on  a particular 
vibrational  level. 

On  Figure  4.21  the  cavity  fluxes  F^  ^ for  a number  of  vibration 
levels  and  the  total  flux,  are  given  for  the  same  conditions.  From  this 
result  the  total  output  power  Is  calculated  to  be  54  kj/lb  with  most  of  this 
power  radiated  from  the  first  three  vibrational  levels.  As  Is  seen  In 
Figure  4.21  the  cavity  was  allowed  to  extend  arbitrarily  to  300  cm  downstream. 
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Figure  4.19.  Maximum  small  signal  gain  for  the  v=l  to  7 vibrational 

levels  of  HCl.  Critical  gain  is  calculated  for  i=100  cm, 
R=0 .98,  and  L=0 . 2 . 


SATURATED  GAIN 


Figure  4.20.  Saturated  gain  for  v=l  to  4 vibrational  levels  of  HCl. 

Cavity  conditions:  Jl=100  cm,  R“0.98,  and  L=0.2. 
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Figure  4.21.  Calculated  inttacavity  flujr  for  v - 1 
to  3 vibrational  transition  of  HCl. 
Also  shown  is  the  total  flux,  i.  e.  , 
the  sum  of  all  levels,  including  those 
not  shown. 
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However,  the  bulk  of  the  energy  is  extracted  early  in  the  flow.  For  example,  at 
X = 100  cm  the  output  power  la  already  35  kj/lb.  It  is  emphasized  that  the 
calculation  presented  in  the  figures  Is  only  illustrative,  and  no  attempt 
has  yet  been  made  to  optimize  these  results. 

This  power  extraction  routine  can  also  be  used  with  the  H.O.  model. 

In  this  situation,  there  is  only  a single  vibrational  level  and  the  calculation 
is  correspondingly  simpler.  The  calculated  output  power  in  the  H.O.  approxima- 
tion for  the  same  assumed  conditions  as  were  shown  previously  for  the  A. H.O. 
case  is  given  in  Figure  4; 22.  A comparison  of  Figure  4.22  with  the  A. H.O.  re- 
sult indicates  that  the  H.O.  calculation  yields  about  twice  the  extracted  power 
as  the  A. H.O.  result.  This  illustrates  one  severe  limitation  of  this  simple 
H.O.  model  of  a vibration  transfer  laser. 


4 . 5 Advanced  H„-HC1  Gas  Dynamic  Laser  System  Characteristics 

The  end  objective  of  the  work  performed  is  to  develop  a high  energy 
laser  device.  The  work  performed  during  this  phase  is  now  sufficient  to  deter- 
mine the  weight  and  volume  characteristics  of  a high  energy  laser  system. 

A typical  HEL  device  size  was  selected.  The  run  time  for  the  device  consists 
of  20-three  second  bursts  for  a total  run  time  of  60  seconds.  The  weight  and 
volume  of  the  laser  system  depends  on  the  propellant  weight  yield  and  the 
specific  power.  The  specific  power  used  to  estimate  the  system  size  is  100  kj/lb. 
This  value  represents  an  upper  limit  which  depends  on  the  effect  of  H-atoms  on 
specific  power.  The  H-atom  problem  will  be  addressed  in  much  greater  detail 
in  the  next  year. 

The  laser  device  has  been  divided  into  three  basic  components  which 
are:  (1)  the  combustor;  (2)  the  nozzle  bank;  and  (3)  the  cavity  .and  diffuser. 

A sketch  of  the  device  is  shown  in  Figure  4.23.  For  the  purpose  of  weight  and 
volume  estimates,  the  combustor  is  assumed  to  be  comprised  of  1"  thick  steel, 
the  nozzle  bank  is  assumed  to  be  solid  steel  and  the  cavity  and  diffuser  are 
assumed  to  consist  of  a 3"  thermal  skin.  The  thermal  skin  is  comprised  of  a 
.25  inch  thick  steel  plate  outer  shell  and  a .25  inch  thick  steel  plate  inner 
shell  as  shown  in  the  inset  of  Figure  4* 23. 
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Figure  ,4.22.  Calculated  output  power  for  H.  O. 

calculation.  Cavity  conditions 
= 100  cnu,  R = 0,  98  and  L = 0.  2 


0.25  in.  Thick  Steel 


0.25  in.  Thick  Steel 

Thermal  Skin 


Figure  4.23.  Laser  Device 
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The  weight  and  volume  of  the  components  is  shown  in  Table  4.5  for 
each  of  the  three  candidate  propellants.  The  overall  envelope  of  the  device 
is  obtained  using  a packing  factor  of  2.  The  weight  and  volume  of  the  indi- 
vidual gas  generators  are  shown  in  Table  4.6.  The  gas  generator  weight  and 
volume  estimates  include  a filter  system  for  propellant  candidates  2 and  3. 

A more  detailed  breakdown  of  the  gas  generators  is  shown  in  Table  4.7. 

The  overall  system  weights  and  volumes  are  shown  in  Table  4.8. 

Also  shown  in  Table  4.8  are  the  specific  weight  power  of  the  laser  devices, 
the  specific  volumetric  power  and  the  specific  nozzle  exit  area  power.  The 
best  device  performance  is  obtained  from  the  HDB  propellant  system,  whereas 
the  other  two  systems  have  about  the  same  system  weight  and  volume.  Depending 
on  whether  the  application  of  the  laser  is  weight  or  volume  limited  either 
the  device  specific  weight  or  volume  parameter  should  be  selected  as  the  basis 
for  comparison  to  other  laser  systems. 
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TABLE  4.5 


WEIGHT  AND  VOLUME  OF  COMPONENTS 


PROP.  1 and 

ITEM 

WEIGHT (lb) 

COMBUSTOR 

330 

NOZZLE 

378 

CAVITY  AND  DIFFUSER 

836 

TOTAL  1,544 


t (L«2  ft) 

PROP.  2 

(L=2.5  ft) 

VOLUME (ft^) 

WEIGHT (lb) 

VOLUME (ft^) 

1.6 

420 

2 

0.73 

500 

1 

12.0 

1,040 

16 

14.33 

1,960  . 

19 

28.66 

38 

PF  = 2 


TABLE  4.6 


WEIGHT  AND  VOLUME  OF  GAS  GENERATOR 


PROPELLANT 

1 

2 

3 


WEIGHT  (lb) 
4,900  (245) 
10,100  (506) 
11,580  (579) 

( ) DENOTES  PER  SHOT 


VOLUME  (ft^) 
45  (2.24) 

88  (4.4) 

104  (5.2) 


132 


PROPELLANT 


(H^+HCl) 

Comp  H„/HC1  Yield 

ID  Composition,  Wt  % Molar  Wt  % Vol  % °K 


1 11.821%  H WO  33.3  10.8  95.9  2808 

33.475%  Zr 
38.704%  HDB 
16%  AP 


2 10.482%  NH.NO  38.2  3.84  89.5  2799 

47.132%  ZrH 
38.186%  Zr(OH) 

4%  AP 


3 40.879%  32.6  2.59  98.0  2799 

42.968%  ZrH 
13.153%  Zr(OH) 

3%  AP 


TABLE  4.7 
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WEIGHT  AND  VOLUME  OF  SYSTEM 


PROPELLANT 

1 

2 

3 


WEIGHT (ft) 
6,444 
12,060 
13,124 


I 


I 


i 


TABLE  4.8 


VOLUME (ft^) 

kj/lb 

MZfti 

k.i  /in^ 

73 

0.062 

5.47 

4.7 

126 

0.033 

3.17 

3.2 

132 

0.030 

3.03 

4.7 

5.0 


CONCLUSIONS 


The  results  of  the  Phase  I study  Indicated  that  an  advanced  H2-HCI 
gas  dynamic  laser  was  indeed  feasible.  The  objective  of  the  current  phase, 
Phase  II,  was  to  experimentally  verify  the  V-V  and  V-T  kinetics  for  the  impor- 
tant H2-HCI  processes  and  to  develop  solid  propellant  gas  generators  that 
could  be  used  to  generate  the  required  gases  at  attractive  weight  yields. 

The  results  of  the  current  phase  have  verified  that  the  kinetics  are  capable 
with  operation  of  a H2-HCI  laser.  In  addition,  the  gas  generator  development 
program  has  succeeded  in  producing  the  required  gases  at  weight  yields  that 
make  the  system  attractive  for  Army  field  application.  The  potential  of  even 
higher  gas  yields  that  would  make  the  system  even  more  attractive  is  also 
possible  but  requires  additional  propellant  development.  The  fundamental 
remaining  problem  area  is  the  effect  of  H-atoms  on  performance.  This  problem 
area  is  currently  under  investigation  and  represents  the  only  remaining 
barrier  that  prevents  the  development  of  an  advanced  H2-HCI  gas  dynamic  laser. 
It  is  believed  that  this  problem  may  have  to  be  resolved  by  the  use  of  an 
adequate  H-atom  scavenging  technique. 
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APPENDIX  A 


Theoretical  values  of  flame  temperature  and  combustion  product  compo- 
sitions are  tabulated  here  for  selected  propellants  as  a function  of  propellant 
compositions.  Values  of  flame  temperature  are  for  adiabatic  combustion  with 
chemical  equilibrium  at  750  psla  (generally).  Values  for  combustion  product 
compositions  are  for  throat  conditions,  assuming  expansion  of  only  the  gaseous 
species  that  exist  at  chamber  conditions.  The  expansion  is  assumed  to  be 
isentropic,  with  chemical  equilibrium  throughout.  The  calculations  are  com- 
puterized results,  using  JANAF  thermochemical  data. 

It  should  be  noted  that  Al202(c)  and  BN(c)  appear  as  combustion 
products  at  throat  conditions  of  all  propellants  containing  A1  (or  AlH^)  and 
HDB,  respectively.  This  phenomenon  has  been  discussed  in  the  text  and  is 
essentially  the  result  of  condensation  reactions  of  gaseous  species  during 
expansion.  The  effect  is  not  large,  amounting  to  *^2%  by  weight  of  combustion 
products  appearing  as  condensed  at  the  throat. 

In  the  tables,  the  gas  composition  is  expressed  in  terms  of  mol/100  gm 
of  propellant.  Some  of  the  tables  show  the  column  entitled  "Eng"  which  refers 
to  the  total  moles  of  gaseous  combustion  products  at  the  throat.  The  molar  ratio 
H2/HCI  is  seen  to  be  quotient  of  the  entry  in  Column  7 (H2)  divided  by  that  in 
Column  8 (HCl).  The  yield  (H2  + HCl)  is  seen  to  be  the  sum  of  the  products 

2.016  (Col.  7)  + 36.461  (Col.  8) 

for  the  expression  for  weight  percent.  The  yield  (H2  + HCl)  for  volume  percent 
is  the  molar  content  of  (H2  + HCl)  in  the  combustion  products , and  is 

(H2  + HCl) /Eng. 


For  most  of  the  tables,  this  quantity  1s  derived. 
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Propellant  System 

, table 


Oxidizer 

Fuel 

Coolant 

Chlorides 

Number 

NH,NO„ 
4 3 

ZrH2 

Zr(OH)^ 

NH,C1 

4 

A1 

NH^CIO^ 

A2 

(CH  ) 
z n 

NH,  Cl 
4 

A3 

NH^CIO^ 

A4 

. 

HTPB 

NH^CIO^ 

A5 

Zr 

HTPB 

NH,C10, 
4 4 

A6 

HDB 

NH^Cl 

A7 

NH^CIO^ 

A8 

HDB/CH2 

A9 

A1 

AKOH)^ 

NH,C1 

4 

AlO 

NH^CIO^ 

All 

(CH  ) 
z n 

NH,C1 

,4 

A12 

NH.CIO,  A13 

4 4 


HTPB 


A14 


AIH3 

A1(0H)3 

NH^Cl 

A1(0H)3/(CH2)^ 

A15 

A15 


N2H^N02  ZrH2  Zr(OH)^  NH^Cl  A16 

N2H^C1  A17 

NH.CIO,  A18 

4 4 

(CH„)  NH,C1  A19 

2 n 4 

N2H^C1  A20 

NH,C10,  A21 

4 4 

HDB  NH,C1  A22 

4 

Zr  ^2^5^^ 

NH,C10,  A24 

4 4 

HDB/Zr(OH)^  ZrCl^  A25 
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Oxidizer  Fuel Coolant Chlorides Number 


' N2H3N03 

A1 

A1(0H)3 

NH^Cl 

A26 

N2H3CI 

A27 

NH^CIO^ 

A28 

<™2>n 

■ NH.Cl 
4 

A29 

N2H3CI 

A30 

NH^CIO^ 

A31 

HDB 

NH^Cl 

A32 

AIH3 

A1 (OH) 3 

NH.Cl 

4 

A33 

(™2>n 

A34 

A1(0H)3 

NH^CIO^ 

A 35 

Zr(OH)^ 

NH^Cl 

A36 

■■■ 

NH.CIO, 
4 4 

A36 

NH^CIO^ 

A1 

HDB 

A37 

Zr 

A38 

A1 

* 

ZrCl^ 

A38 

Zr 

Zr(OH)^/(CH2)j^ 

A39 

ZrH2/Zr(OH)^ 

A40 

H2W0^ 

ZrH2 

Zr(OH)^ 

NH4CIO4 

A41 

A42 

ZrCl^ 

A43 

ZrCl^/NH^ClO^ 

A44 

Zr 

HDB 

NH^CIO^ 

A45 

HDB/ (CH2)^ 

A46 

HDB 

NH^Cl 

A47 

• 

Zr/ZrH2 

Zr(OH)^/(CH2>^ 

NH.Cl 

4 

' \ 

A48 

NH.CIO, 
4 4 


Zr/ZrH. 


zrn^lzxion)^ 


N2H^C1 

NH.CIO, 
4 4 


NH,C1 

4 

NH,C10, 
4 4 


Zr0(N02)2-2H20  ZrH2 


Zr(OH), 


(CH2)„ 


NH.Cl 

4 

NH,  CIO, 
4 4 

NH^Cl 

NH,C10, 
4 4 


Hi 


